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►^ ■ Abstract: It has been recently demonstrated that the thermal partition function of any 

1/^ ■ large A^ Chern-Simons gauge theories on 5^, coupled to fundamental matter, reduces to a 

(N 
!>■ , 

^T) , Chern-Simons theories at temperature T were determined in ||I|. The large N saddle point 

Tjlj- I equations for these theories were determined in the same paper, and were solved in the low 



temperature phase. In this paper we find exact solutions for these saddle point equations 
in three other phases of these theories and thereby explicitly determine the free energy of 
the corresponding theories at all values of T^/A. As anticipated on general grounds in Q, 
our results are in perfect agreement with conjectured level rank type bosonization dualities 



Vh ■ between pairs of such theories. 
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1. Introduction 

Recently, three dimensional field theories have been fascinating in the context of AdS/CFT 
(or dS/CFT) correspondence which do not necessarily rely on the supersymmetry H. 
It has been recently conjectured that level k Chern-Simons theories with matter in the 
fundamental/bifundamental representations admit a dual description governed by parity 
violating Vasiliev's higher spin equations |^ ^ (See P for more detail and references) at 
finite values of the 't Hooft coupling A = -^ [^i 01 • 

Sometimes, analysis of boundary field theories could be applied to study bulk grav- 
ity theories through the AdS/CFT correspondence when the direct analysis of the gravity 
theory is difficult. For example, it is known that the AdS/CFT correspondence maps 
deconfinement transitions of large A^ gauge theories on spheres to gravitational phase tran- 
sitions involving black hole nucleation |^]. This observation has motivated the intensive 
study of the phase structures as well as the deconfinement phase transitions of large A^ 
p-|-l dimensional Yang Mills theories (coupled to adjoint and fundamental matter) on S^. 
Now, we can expect similar application of the Chern-Simons matter theories to know the 
phase structures of parity violating Vasiliev's higher spin gravity theories. 

From this motivation, in |||], it was studied the finite temperature phase structure 
of renormalized level k U(A^) Chern-Simons (CS) theories coupled to a finite number of 
fundamental fields on S"^ in the t' Hooft limit A^ — )• oo, /c — )• oo with A = -^ fixed. The phase 
structure was studied by using a simple toy large A^ Gross- Witten-Wadia (GWW) matrix 
integral, and it turned out there is a rich phase structure caused by the summation over the 
U{1) flux sectors in the Chern-Simons theory. This summation resulted that the eigenvalue 
must be discretized and then the eigenvalue density function must be saturated from upper 
bound ^^. The saturation is first suggested by Q. This discretization generates the new 
phases so called the "upper gap phase" and the "two gap phase" which were absent in 
Yang- Mills theories. ^ And the level-rank duality has been confirmed in the GWW toy 



model. (See |14| , 15 1 for a relatively recent discussion of level-rank duality and references to 
earlier work.) But the phase structures of the actual Chern-Simons matter theories have 
not been worked out. Moreover we have not had the perfect analytic proof of the duality 
even in the GWW model. 

Then in this paper, we will work out the study of the phase structure of several Chern- 
Simons matter theories: (1) the CS theory minimally coupled to the fundamental fermion 



in, (2) the CS theory coupled to massless critical bosons [^, 17], and (3) the J\f = 2 
super symmetric CS theory with a single fundamental chiral multiplet ||l8| , |l9| ]. (SUSY CS 
matter theory) We will also provide the perfect analytic proof of the level-rank duality 
(Giveon-Kutasov type duality) on the 5^ x S^ between the regular fermion theory and the 
critical boson theory, and the level-rank self-duality in the SUSY CS matter theory as well 
as the GWW matrix integration. In the proof, we supply several formulae which can be 



""^The same upper bound for the density function appears in two dimensional Yang-Mills theory {p = 1), in 
which case the situation becomes similar due to the fact that there is no propagating degrees of freedom for 
the gauge field. A phase transition relevant to this upper bound was studied in 2d (q-deformed) Yang-Mills 
theory on S^ tld, O, [12, [13|, which we will see from Chern-Simons theory. 



used for any related analysis to the duality in general CS matter theories on S'^ x S^. 

1.1 Outline of this paper 

The rest of this paper is organized as follows: In section 2, as a preliminary, we will 
review the Landau-Ginzburg description of the partition function of Chern-Simons theories 
coupled to fundamental matters. And we will supply techniques to evaluate eigenvalue 
densities as well as partition functions by using complex line integration in the complex 
plane. In section 3, we will elaborate the phase structure of the regular fermion theory 
and the critical boson theory. And then we examine the level-rank duality between them 
and we will confirm the agreement. In section 4, we will study the phase structure of the 
M = 2 super symmetric CS matter theory and examine the Giveon-Kutasov type duality. 
We will confirm the self-duality of the theory. In section 5, we will give the analytic proof 
of the level-rank self-duality in the two gap phase of the GWW toy model. Section 6 is 
devoted to the summary and discussion. 

The content of the appendix is as follows: We put several definitions of the cuts region 
and cut functions in appendix ^. And we give the proof of important formulae used in 
the proof of the duality and others in appendix g In appendix ^, we provide the detailed 
analysis of the large C behavior of the eigenvalue density functions. In the appendix ^, 
we have also proved that the eigenvalue distributions at ^ — )• oo converge to the universal 
distribution ( |2.16| ). 

2. Preliminaries 

2.1 High temperature effective action of the Chern-Simons-fundamental matter 
theory on S"^ x S^ 

We will consider the level k U{N) fundamental matter Chern-Simons theories in the 't 
Hooft large A^ limit, with^ 

V2T^ = NC, \ = j (2.1) 

where C and A are held fixed in the large N limit. The theories have recently studied 
intensively at finite A [^, 16, |^, 23, ^, |^, ^, |2^, 26, 27, 28, ^, and thermal partition 



functions of the theories are studied in ||2^, p, |7| . 

We will consider the partition function of the theories 



Zcs= f DADfi e'^'^'^i^'^^+l^')-^-''^^ (2.2) 



where Dfi denote the integration measure of matter fields. To calculate it, we will integrate 
the matter fields first, and obtain the effective action which depends on the gauge fields. 



^We use the dimensional reduction regulation scheme through this paper. In this case |fc| = |fcyA/| + TV, 
where kvM is the level of the Chern-Simons theory regulated by including an infinitesimal Yang Mills term 
in the action, and k is the level of the theory regulated in the dimensional reduction scheme fed, fell . For 
this reason this class of Chern-Simons theories may be well-defined only in the range jA| < 1. Hereafter we 
assume the 't Hooft coupling is always positive for simplicity. Our results are easily generalized to the case 
A is negative by taking the absolute value of it. 



According to the study in |1|, the partition function is given as 

Zcs = /"l)^e^^Tr/(AdA+§A3)-5.,,(^) (2.3) 

with Seff{U) which is a function of the two dimensional holonomy fields U{x) around 
the thermal circle 5^. Here x of U{x) indicates a point on the 5^, and its eigenvalue is 
described as e*"'"^^) where m is the index of N x N matrix. The effective action Seff{U) 
would be obtained by summing up all vacuum graphs in which involve at least one matter 
field. The graphs which do not involve any matter fields are not summed yet, but they are 
summed at the stage of the path integration in (2.3). 

To calculate (^^), it is useful to observe that Seff{U) is ultralocal in large A^ as well 
as in high temperature. The Seff{U) would be expanded by a series of local operators, 

SeffiU) = I (fx (T^^viU) + vi{U)TrD,UD'U + ...). (2.4) 

Please note that the scaling (^^ converts higher temperature expansion (^^) into expan- 
sion in inverse power of N. Here the first term is leading 0{N'^) and the second term is 
0{N^) and terms . . . are further suppressed at large N. ^ Then at large N, the effective 
action is simplified to be the leading term 

SeffiU) = jd^xT^^viU). (2.5) 



This (2^) is originally suggested in [^ , which is a consequence of the entire effect of matter 
loops on gauge dynamics at temperature yN and at leading order in N. By ( p.SD , the 
partition function would be given as 

^l^^-T^Jd^x^viUix))y^ (2.6) 

= (e-^^^^^(^));v,. = (e-^^^(^));v,. 

where 

(^)iv,fc (2.7) 

is the expectation value of ^ in the pure U{N) Chern-Simons theory at level k. Since the 
Chern-Simons theory is topological, expectation values become independent of x, then the 
X dependence is removed in the step from the second line to the third line of ( |2.6D . Here 
V2 is the volume of S"^. 

2.2 Path integration to obtain the partition function, deconfinement phase 
transition and further phase transition caused by flux 



We can perform the calculation of (|2.6| ) in the same manner as 1 3C ] . You can see the detail 



of the way how to compute the path integration in m . With keeping the effective potential 



^Every term in the effective action (2.4) is of order A'^ at fixed T, as tlie action is generated by integrating 
out fundamental fields. 



term N(^v{U) intact, we achieve the foUowing value 



Zcs= ff[daj m 2sin 



am{nm) -ai{ni) \ \ ^-ncv{u) 



E 



^ikMjOj 



(2.8) 



Oim{nm) - ai{ni) 



-NCviU) 



y 6{kaj — 2'Kn) 



nez 



where Mj is a constant units of flux in the U{l)j (j = 1, . . . N) factor. As we can see in 

Iu€ 

(2. 



the last line, the summation over flux Mj in e*'^^^^"^ constrains aj to take discrete values 



27rn, 



Oii 



Then the partition function would be 

N oo 



zcs=n E 



Y\ 2 sin 



ai{n) - amin] 



'NCo(U) 



(2.10) 



where the summation over rim is restricted so that no two nm are allowed to be equal. 

2.3 Eigenvalue density 

We will estimate ( p. 10 ) in the 't Hooft limit, k — )• oo, A^ — )• oo, A = -^ : fixed. The summa- 



tion over the eigenvalue ( 2.10| ) in the limit is dominated by the saddle point configuration 
of the eigenvalues at large N. The saddle points minimize the following potential 



y([/)-^ln2sin 



am - Oil 



where V{U) = ^ t;(a„) = ^ e*° 



(2.11) 



where V{U) = NC,v{U). The saddle points are obtained by solutions of the following 

equation 

otm - ai 



V'{am) = y^ cot — 

In terms of eigenvalue density function p{a), the equation would be represented as 

V'{ao) = NV [ da cot "°~" p(a) 



(2.12) 



(2.13) 



where p{a) = ^ YlZ=i '^(« - «m)- 

We can see that this saddle point equation problem is very similar to the Gross- Witten- 
Wadia (GWW) problem |3l|, |3|, |3| in the Yang-Mills theory on SPxS^. In a usual GWW 
problem in Yang-Mills theories, a thermal partition function obtained after integrating out 
all massive modes is given by the integration over the single unitary matrix as 



N 



Zym = J_| 



m=l 



dar 






sm|^i^]e-^--(^) 



(2.14) 



and there is a competition between potential VVm(C^) that tends to clump eigenvalues, 
and the measure factor (Vandermonde determinant) which tends to repel them. ^ In the 
low temperature with small C, the repulsive force by the measure factor is stronger. Then 
the eigenvalue density would have support everywhere on — vr < \a\ < vr. But on the other 
hand, in the high enough temperature, the attracting force caused by the potential Vym{U) 
becomes stronger than the repulsive force, and the eigenvalue will be clumped. Then in 
the high temperature the eigenvalue distribution has support only on a finite arc, there 
would be a domain of eigenvalue so called "lower gap" such that {a\a ~ a + 2tt, p{a) = 0}, 
and phase transition would occur. And in the extremely high temperature T — )• oo the 
eigenvalue density function would be clumped to be a delta function p{a) ~ 6{a). 

Also in the current Chern-Simons case ( ^.lOl ), such a competition between V(U) and 



measure term exists. But in this case, different from the usual Yang-Mills case ( 2.14] ), there 



is a constraints on the eigenvalue ( |2.9D . ( p. 9]) saturates the eigenvalue density from above 
as 

while the eigenvalue density in the Yang-Mills theory is not bounded from above (is bounded 
only from below). Saddle points of ( 2.14| ) do not always obey the inequality ( 2.15| ), then the 



Yang-Mills solutions violating ( p. 151) will not be the solution of the Chern-Simons saddle 



point equations ( p. 13 ) and ( |2.10| ). Instead, the current Chern-Simons theory admits new 



classes of solutions saturating the upper bound of the inequality over several arcs along the 
unit circle. Hence in this Chern-Simons case, there can exist not only the "lower gaps" but 
also the "upper gaps", which are arcs over which the upper bound of (|2.15|) are saturated. 



Hence the phase structures of the current Chern-Simons theories ( 2.101 ) would be dif 



ferent from the one of usual Yang-Mills theory (2.14). In a usual Yang-Mills theory case, 
there are only two phases, "no gap phase" and "lower gap phase". In the no gap phase, the 
eigenvalue has support everywhere on the unit circle in the complex plane, on the other 
hand in the lower gap phase, the eigenvalue has support only on a finite arcs. On the 
other hand, in the current Chern-Simons case, due to the existence of upper bound ( 2.15| ), 



not only the no gap and lower gap phases, but also "upper gap" and "two gap" phases 
exist. In the upper gap phase, there is one upper gap where the upper bound of (l2.15| ) is 
saturated, and the eigenvalue density has support everywhere on the unit circle. In the two 
gap phase, there are one upper gap but also there is one lower gap where the eigenvalue 
density vanishes. 

In the Chern-Simons case, the eigenvalue distribution is in the no gap phase in the 
low temperature, and if we increase the temperature, eigenvalue density starts to clump 
to a = 0, and then we can expect lower gaps or upper gaps will show up. In the small 't 



Hooft coupling region A < Ac, since the upper bound in (2.15) is large enough, then lower 
gaps will show up before the maximum of the eigenvalue density function achieves to the 
upper bound. Then it transits to lower gap phase first. And if we increase the temperature 



*The efTective potential Vym{U) was computed in free gauge theories |34l, B5|; it has also been evaluated 
at higher orders in perturbation theory in special examples |3q, p7^, pq, p9|. At least in perturbation theory 
Ipq] and perhaps beyond Wd, kl|, the potential Vym{U) is an analytic function of U. 
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Figure 1: The eigenvalue distribution p{a) in the no gap phase Fig. 1(a), in the lower gap 
phase Fig. |l(b)| , in the upper gap phase Fig. 1(c), and in the two gap phase Fig. |l(d)| . These are 
the graph calculated in the GWW toy model. We can see that the upper gap one Fig, 1(c) at 
0.51 is saturated from above as p{a) = 2xSbiT: ^* ^^^ upper gap region \a\ < 0.4. We can also 



A 



see that the two gap density function Fig. 1(d) has both a lower gap and an upper gap 



further, the maxiinuni of the eigenvalue density p eventually reaches to the upper bound 
■^:^ and it transits to the two gap phase. On the other hand, in large A with A > Ac, the 



upper bound in (2.15) is small, then the maximum of the eigenvalue reaches to the upper 
bound before the lower gap shows up. Then in the large A region, it transits to upper 
gap phase first. If we increase the temperature further, since the eigenvalue density keeps 
clumping, lower gaps will show up and it transits to two gap phase. In the large C limits, 
the eigenvalue density will be universal configuration 



p{a) 



1 



2nX ^"""<"^) 
(|a|>7rA) 



(2.16) 



which is the nearest thing to a 5 function permitted by the effective Fermi statistics of the 
eigenvalues < p{a) < 2:^- This is in perfect agreement with the results of |p in which 
they used the Hamiltonian methods J42D . 

For your reference, you can see the graph of the eigenvalue density function in each 



phase listed in Figs. 1(a) ~ 1(d) , And you can also see the phase diagrams of Chern-Simons 



matter theories in Figs. 2(a) 4(a) and 5(a) 



2.4 How to calculate the eigenvalue density in lower gap, upper gap and two 
gap phases 

Here we will review how to obtain the eigenvalue density in the lower gap, the upper 
gap and the two gap phases. To calculate the eigenvalue density, we will decompose the 
eigenvalue density to p{a) = po{a) + ip{a). And substitute it to the saddle point equation 
and rewrite it as 

NV I da V'(a)cot ( "°~" ) = U{ao) 



da V'(o) = A[po] 
U{a) = V'{a) -NV de po{9) cot " 



A[po] = 1 - da po{a). 



Here "0(0^) is nonzero in the compliment of both upper gaps and lower gaps. Starting from 
this, we turn to the complex variables 

Note that 

dz uq — a .zq + z 

da = — , cot = I . 

iz 2 zo — z 

Here ( 2.17]) will be rewritten in terms of these variables as 



,r^ f dz zo + z ^ , . ^^, , 
NV / ^ ^{z) = U{zo), 

J Z Zq- Z 

dz , , . . 
— ip{z) = A 
iz 



(2.18) 



where the integrals in ( |2.18| ) run counterclockwise over the unit circle in the complex plane. 



Suppose that the solution to this equation, '<p{z), has support on n connected arcs on the 
unit circle in the complex plane. We denote the beginning and endpoints of these arcs by 
Ai = e*"* and Bi = e*^' (i = 1 . . . n). These are the endpoints of lower gaps or upper gaps 
as well. Our convention is that the points A\, Bi, A2, -B2 • • • An, Bn sequentially follow 
each other counterclockwise on the unit circle. We refer to the n arcs {Ai,Bi) as 'cuts'. 
The arcs {Bi, Aj+i) (as also {Bn, Ai) ) are referred to as gaps. By definition, ^ has support 
only along the cuts on the unit circle. 

We use the (as yet unknown) function 'ip{z) to define an analytic function ^{u) on the 
complex plane 

where the integral is taken counterclockwise along the n cuts, {Ai, Bi) of the unit circle in 
the complex plane. It is easy to see that the analytic function <l>(u) is discontinuous along 



these n cuts. Let ^{z)~^ denote the hmit of $(u) as it approaches a cut from \u\ > 1, and 
let ^~{z) denote the hmit of <5(n) as it approaches a cut from |n| < 1. Then ^ 

$+(z) - $"(z) = 47rV'(z) (2.20) 

and 



<!>+{z) + <^-{z) = 2^V -. p(uj) = —^ (2.21) 

^-^ J A, lUJ Z- UJ iN 



"^^ dujz + uj ^_^_ 2U{z) 
I Ai iu) Z — OJ 

(the first equahty is the definition of the principal value, while the second equality follows 
using ( 2.1^ )). Moreover it follows immediately from ( p. 18 ) and ( p. 19 ) that 



lim $(ti) = A. (2.22) 

We are now posed with the problem of determining $(«) given its principal value along 
a cut. To determine the ^{u), we introduce following "cut function" h{u) as 

h{u) = V(^i - u){Bi - u){A2 - u){B2 -u)...{An- u){Bn - u). (2.23) 

We define h{u) to have cuts precisely on the n arcs on the unit circle that extend from Ai 
to Bi. This definition fixes the function h{u) up to an overall sign. This sign will cancel 
out in our solution for $ below, and so is uninteresting. For future use we note that when 



h{u)=YY^e 2 e (sm ( ^ 1 - sm (- ^ 11. (2.24) 

m=l ^ ^ ^ ^ ^ ^ 

We use the function h{z) to define a new function, H{z), via the equation 

$(z) = h{z)H{z). 
Using the fact that h^{z) = —h~{z) along the cut, (|2.2lD turns into 

And here we assume that H{u) is holomorphic except the cut region. From ( |2.22| ) and 



d^H), H{v) = 0{^) at large v so that 



Hiv) 

dv / ^ , = (2.26) 

Coo 27r<f - u) 

where the contour Coo runs counterclockwise over a very large circle at infinity. Since H[v) 
is holomorphic except the cut, from the Cauchy's theorem, 

0=/ dv^P^-f dv^P^-Hiu) = -{ dv^P^-H{u). 
r 2'Kiiv — u) fr , 2'Kiiv — u) r t 2TTi{v — u) 

(2.27) 



^We use analogous notation for other analytic functions below. Completely explicitly, let z denote a 
complex number of unit norm. The symbol F^{z) will denote the limit of F{u) as u -^ z from above (i.e. 
from \u\ > 1), while F~{z) is the limit of the same function as u — > 2 from below (i.e. from |m| < 1). Note 
that along a cut F~ {z) — —F'^{z). 



10 



The contour Cents are loops enclosing each of the n cuts {Ai,Bi), but the loops do not 
enclose the point u. Then from ( 2.25| ), 



H(u) = -d. d. ^(-) .=1/ ...... ^(f) .=J-/ d.- ^(^) 



/a„.. 27ri(t;-tx) 7r7i_ Nh+{z){z-u) 2t, Jc.^^^ Nh{v){v - u)' 

(2.28) 
In the equation above, f is a variable on the complex plane while 2; is a variable on the 
unit circle of the complex plane. The integration region Lares runs counterclockwise along 
n cuts on the unit circle (This is not loop !), integration region for each cut is from Ai to 
Bi. The third equality uses h^{z) = —h~{z) together with the assumption that U{z) has 
no singularities on the n cuts. 

Now by applying the Cauchy's theorem for the integrand ^y^. wj_^v again, we will 
obtain the following 

H{u) = i- / dv ^(") 



2vr/a„,. Nh{v){v-u) 

2^fcJ NKz){z-u) Nh{u) 2^^''^''^=^-Nh{z){z-uy 

The last term is the sum of the residues of U{z). Here we have assumed that U{z) is 
meromorphic function of z. 

Based on H{u) in ( p. 29 ), and from $"'"(tt) — $~(u) = 47r'0(ii), by combining with 



p{a) = il^{a) + po{a)^ we can obtain the eigenvalue density p{a). From next section, we 
will obtain the phase structure of the regular fermion theory, the critical boson theory 
and the M = 2 supersymmetric CS matter theory by using the techniques supplied in this 
subsection. 

3. Higher temperature phases of the regular fermion theory and critical 
boson theory and the duality 

3.1 Regular fermion theory 

In this subsection we study the level k U{N) Chern-Simons theory coupled to massless 
fundamental fermions. The Lagrangian of the theory is presented in equation (2.1) of Q. 
From the equation (3.5) of [|^, effective potential V{U) is obtained as 

V{U) = —^(j-e' + 3 dap{a) I dy y(ln(l + e'^^"^") + ln(l + 6"^+*")) J 

^V'-f[p,N-c,C], 

(3.1) 

where c determines the thermal mass of the fermions.^ The value of c is obtained by 
extremizing V{U) w.r.t c at fixed p, (, i.e. c obeys the equation 

/I c ~\~ tot c — tiX \ 

da p{a) ( In 2 cosh(^ — ) + In 2 cosh( — - — ) J . (3.3) 



'Et ~ (?T^ is the thermal mass of the fundamental fermions. More precisely, the fermionic self energy 



11 



c 


\ 










\ 

Lower gap 


Two 


gap 






Phase 


\ 




/ 




{x,i;) = (KQ~~~~~~ 


\ 


■<^ 


/ 

Upper 
^ap 




No gap 







(a) 

Figure 2: Phase diagram of the regular fermion theory. Here (Ac, Cc) = (0.596967, 2.86454) is the 
quadruple phase transition point where the four phases (no gap, lower gap, upper gap, two gap 
phases) coexist. 



The general form of the free energy of the regular fermion theory on S*^ x S^ is given as 






V''-^{p,N\-N^V 



da 



df3 pia)p{f3) log 



2 sin 



a — 13 



--V'-f[p,N] + F2[p,N]. 



(3.4) 



To obtain the free energy in each phase, we only have to evaluate the eigenvalue density p 
and c in each phase and just substitute into (3^). Please note that c can be determined if 
the eigenvalue density is determined. So obtaining the eigenvalue density in each phase is 
equivalent to obtaining the free energy in each phase. 

For later use, we will give the form of V'{z) here. We obtain V'{z) from ( |3.lD as 



V'{z) 



2lT 



dy y 



-le 



z + e-y 



+ 



le 



+ e-y 



(3.5) 



The phase structure of this theory is depicted in Fig. 2(a), as obtained in following 
subsubsections. 

3.1.1 Lower gap phase 

To obtain the eigenvalue density in the lower gap phase, we will employ the cut region and 
the cut function described in appendix |A.l . The function H{u) as well as <I>(n) = h{u)H{u) 
in the lower gap phase are obtained by using ( |2.29 ) with substituting pQ{a) = 0. Because 
/>o(a) = 0, U{z) in (|2!29| ) is identical to V {z) as 



U{z) = V'{z) 



NC 
2tt 



dy y 



-le 



z + e-y 



+ 



le 



is given by Et(p) = /(/3ps)pa/ + ip giPps)^ , where 

da p{a) I In 2 cosh( 
V 

9iy) = ^~ fiyf- 



n,)-'- 



Vp + fi^ + i 



+ In 2 cosh( 



+ e-y 



Vy^T^ 



(3.6) 



(3.2) 
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By substituting the above into ( 2.29| ), we immediately obtain H{u) as well as ^{i 
h{u)H{u) as, 






yh{u){u + l){l + eV) 






+ ir{ '^y y — r~^ ~ / ^y y -i , -u 

2tt \J^ u + e y Jc u ^ + e y 

From $+(n) — $~(n) = 4:iTp{u) at the unit circle u = e*" with —b < a < b, we obtain 
the eigenvalue density in the lower gap phase as 



C /„■ 2 ^ • 2 '^ /"°° J y cos ^ cosh I 



p{a) = I- A /sin - — sin — ■ / dy 



v27r2 V 2 2 Jg (cosh y + cos a) y^ (cosh y + cos 6) (3.8) 



And at vr > |a| > 6, /9(a) = 0. By substituting the eigenvalue density ( |3.8| ) into ( p^ ) 
and ( |3.4D , we can obtain c as well as the free energy in the lower gap phase of the regular 
fermion theory. 

The condition lim^-^oo ^("u) = 1 requires the following condition 

Mr/(C,c,6) = ^ r dx fi^ - (^ + "^ ^ ^°g" ") = 1. (3.9) 



By ( |3.9| ) and (|3.3p, we can obtain (6, c) as functions of (A,C) as (6, c) = (6(A, C), c(A, C))- 
(|3.8|), (^^) and (3.3) provide a complete set of the solutions in the lower gap phase of the 



regular fermion theory. 



3.1.2 Upper gap phase 



Next we will search for a solution with no lower gap and one upper gap. The domain of 
cut and the cut function in this upper gap phase are defined in appendix |A.2| . To obtain 
the eigenvalue density based on (|2.29|) , we will take /Oo(a) as non-zero, and po(a) = 2:^ at 



the upper gap region —a < a < a. Then by using ( 2.29 ) with (|3.5| ), we obtain H(u) as 
well as <I>(n) as 

$(m) =$;-/(c, a, C; u) + $2° (A, o; u), where 

^'■f. (c, a, C; u) =^- / dy =^ 

^ 27r7g Y^(e2' + e*")(e?^ + e-")(e?' + u){e-y + u) 

C /""^ e~y C /""^ e-S' (3-10) 

+ ir ^y y — r~^ ~ ^ / ^2/ 2/ _i , ,^ , 

27r jg u + e y 2,7: J c u ^ + e y 

$£« (A, a; u) =-^ f du ——}-—-h{u) + iV I du p^{u) "" + '^ 



Uj{u — bj) ' 



where L^gg runs counterclockwise over upper gap region, and Ldr runs counterclockwise 
over unit circle. 
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From ( ^.10 ), and by taking $+(«) — <!> (u) = 47r(p(u) — po(^)) at the cut region u = e* 
with vr > |q| > a, we can obtain the eigenvalue density function as 



/3(a) 



c 



a 



^ , / sm sm 

27rA V2^2\/ 2 2 



dy 



y\ sin ^1 sinh | 



Vcosh y + cos a(cos a + cosh y) (3-11) 



P«/(C,A;c,o;a). 



And in the region \a\ < \a\, p{a) = 2^- To derive (3.11), we have used a formula 

1 ... 1 



2^2A 



duj 



Lugs 



h{ijj){ijj — u 



■h{u) 



27rA 



(3.12) 



where the u = e*" is located at the cut region. This is proved in appendix B.1.2. By 



substituting the eigenvalue density (3.11) into ( p.3[ ) and ( p. 4] ), we can obtain c as well as 
the free energy in the upper gap phase of the regular fermion theory. 

From the condition, lim^^oo ^{u) = 1 — J^^ dapo{a), we obtain the equation 



K/(C,c,a) = |^^^ 



dx 



logx (1 — x) 



logx 



\/x'^ + 2x cos a + 1 



1 



Here we have used the formula 






Lugs iT'K^) 



1 

A 



(3.13) 



(3.14) 



which is proved in appendix |B.1.1| . By ( |3.13D and ( |3.3D , we can obtain {a,c) as functions 
of (A,C) as {a,c) = (a(A, C), c(A, C))- ( [3.11 ), ( [3.13D and (^) provide a complete set of 
solution in the upper gap phase of the regular fermion theory. 

3.1.3 Two gap phase 

Now we will search for a solution with one lower gap and one upper gap. The detail of our 



two cuts and the cut function h{u) is described in appendix |A. 3 

The function H{u) as well as <I>(n) = h{u)H{u) are obtained by using ( p. 2 

$(7/) =</(C, a, 6, c; tx) + ci>f«(A, a, 6; n) 



^t/(C,a,6,c;n) 



'tg 



$^«(A,a,6;n)^— 



dy 



duj 



ye y{u ^ 



Lug 



where 
yh{u)ey 2u + ey + e'^ 
{u + e-y){u-'^ + e-y) h{-ey) {u + ey){u + e'V) 
hiu) f u + uj 

-—— r + tV / du}po{uj)—, r. 

h{uj){uj-u) Jl^^^ u){u-uj) 



(3.15) 



From lim^^oo ^{u) = 1 — J_ da po{a), we obtain the following two conditions, 

1 ( f 

— — T(a,6) =— y--f(a,6,c), where 

47rA 2ir 



dO 



T{a,b)^ r 

J —a 
fca 

y'-f{a,b,c)= dy 



1 



,2 a 
2 



,2 b 
2 



(3.16) 



y^(cosh y + cos a) (cosh y + cos b) 
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1 C . f 

-A{a,b) =1 — -—Q'''''{a,b,c), where 



47rA 47r' 

A. u^ r ^n COS0 

A(a, 6) = / d9' 

' ,2 a „;„2 e /„;„2 b „;„2 e ' (3-17) 



/>oo 

g''-^{a,b,c) = dy y 



2 "^^^^ 2 V 2 '^^^^ 2 



Y^(cosh y + cos a) (cosh y + cos h) 



By using ( |3.16D , ( p.l7| ) and (^I^), (a, 6, c) are determined as functions of (A, Q) as (a, 6, c) = 
(a(A,C),KA,C),c(A,C)). 

From ( ^.151) and by taking <I>+(n) — <I>^(ti) = 'iT:{p{u) — pq{u)) at the cuts \a\ < \a\ < \b\, 
we obtain the eigenvalue density as 

/o(«) =Ptg{oi) = PiigiC, a, b, c; a) + p2,tg(A, a, b; a), where 
p,,,(C, a, 6, c; a) =^^(a, 6; a) j^ ^V ,^^^,^^.y^ (. cosa + coshy 



P2,tgiX,a,b;a) = ^ T{a,b;a) h{a,b,a), 



-F(a, 6, a) eeW (sm - - sm^ -)(sm^ - - sm^ -), 
fr.fia, b; y) =\/{l + 26-?^ cos a + e-2y)(l + 2e-y cos b + e-'^v), 
Ii[a,b; a) 



(3.18) 



'^ (cos 6* - cosa)J{sm'^ f - sin^ f) (sin^ | - sin^ f) 

And p(a) = at 6 < |a| < vr, and p{a) = ^j at < \a\ < a. During the calculation of 
( |3.18 ), we have used ( ^.16 ) also. Here p2,tg has the same functional form as the eigenvalue 



density in the GWW type matrix integration (7.6) in Q. p^ '{ can be regarded as an 
additional term depending on the detail of the theory. By substituing the eigenvalue 
density ( |3.18 ) into ( |3.3D and ( |3.4D , we can obtain c as well as the free energy in the two 
gap phase. Here by taking uj = e*^, T(a, 6), A(a, 6) and \ sin a\T{a,b, a) Ii{a,b, a) can be 
represented by the complex line integral as 

T(a,6) =-Aij '^Tf^^ 

UJ 



^<"'''* = -''*A.,/"SH' <=*■"" 



sina|J^(a, 6, a)/i(a, 6;a) =i \ 



, h+{u) ( 2 1 



h(uj) \{uj — u) u 



These are useful to discuss the level-rank duality later. 

The combination of ( |3.16| ), ( |3.17| ), ( |3.3D and ( 3.18 ) provides a complete set of solutions 
in the two gap phase. 

At the large C limit, as we proved in appendix |C.1.1| , the eigenvalue density approaches 
the universal distribution ( 2.16| ) because c is kept as finite positive quantity at the limit. 
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In the limit, a,b and the eigenvalue density behave as 

a=7rA--, 6 = 7rA + -, 

^ . , , , / sin(7rA) , , /"°° , ?/ . 

e=8sm(vrA)exp^ ^AC)^ ^^ ,,,H, + cosvrA j + ' ' ' ' (3-20) 

1 _i a — a 



p(a)=— cos 1^ 



At appendix C.2.1, we have demonstrated the above behavior ( |3.20| ). As we can see in the 
appendix, the eigenvalue density function is dominated by p2.tg at the large C limit. And 
the form of the eigenvalue density approaches cos^^ \/oi function same form as (7.11) in 
||T| which is the large C, limit of the one of the GWW model. We can also see that the range 
of the domain of the cut e is smaller than the one in the GWW type matrix integration at 
the same value of (A, C)- 

3.1.4 Phase transition points 

Phase transition points between no gap and lower gap Let us consider the behav- 
ior of the lower gap solutions ( |3.8D and (|3.9| ) at the point b = ir which would correspond to 
the phase transition points from the lower gap to the no gap phase. If we substitute b = n 
into (^^), it becomes 

C /"^ , 2 log X C sr^ fl + cn 



This is exactly same as the condition for the phase transition from the no gap to the lower 
gap discussed in Q, which is obtained by substituting a = vr into the third line of (6.33) 
and by requiring /^(vr) = 0. Based on this, if we substitute 6 = vr into (|3.8|), we can see 



C /. 2^ ■ 2« Z""" , y cos f cosh I 



v27r2 V 2 2 Jg (cosh y + cos a) \J (cosh y + cos b) 

^ J^ "^^ ^ V 2(1 - e~y) ~ 4(1 + 6-?'+*") ~ 4(1 + e-y-^^) ) ^^'^^^ 



oo 



1 C v^/ iNn (l + cn 



n=l 



In the last step we have used ( |3.21 ) and the Taylor expansion and 



'^ dyye-y={^-±P\e--^. (3.23) 



The last line of ( 3.22 ) is exactly same as the eigenvalue density function in the no gap 
phase described in the third line of (6.33) in Q. So, at the phase transition points, the 
lower gap solutions are smoothly connected to the no gap phase solutions. 
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Phase transition points between no gap and upper gap Let us consider the be- 
havior of (3.11) and ( f3.13 ) at a = which would correspond to the phase transition points 
from the upper gap to the no gap phase. If we substitute a = into ( 3.13| ), it becomes 



A 27r 



dx 



2 log X C 



1 + x 



^ >:(-!) 



1 + en 



TT 



n=l 



n^ 



(3.24) 



This is exactly same as the condition for the phase transition from the no gap to the upper 
gap phase discussed in Q, which is obtained by substituting a = into the third line of 
(6.33) and by requiring p{0) = 2^- 

Based on this, if we substitute a = into ( 3.11| ), we can see 

pi") = TT— r T^^a/ (sin^ - - sin^ -) / dy 



2ttX 

1 
2^ 

1 

2^" 



V2' 



vr^ 



y sin ^ sinh | 



+ ^ /_ dyy 



27r2 



/ c 

oo 



1 



-y/cosh y + cos a(cos a + cosh y) 
1 1 



a=0 



^^(— 1)" cosna 



2(1 + 6-2') 4(l + e-^/^ 
1 + en' 



A{l + e~y- 



n=l 



n^ 



(3.25) 
In the last step we have used ( p. 241 ) and the Taylor expansion and the formula ( |3.23| ). 



The last line of ( 3.25 ) is exactly same as the eigenvalue density function in the no gap 
phase, described in the third line of (6.33) in |||]. So the upper gap solutions are smoothly 
connected to the ones in the no gap phase at the phase transition points. 

Phase transition from the lower gap to the t^vo gap In the lower gap phase at 
fixed A, the phase transition from the lower gap phase to the two gap phase will occur 
when the maximum of the eigenvalue density /9(0) reaches to /9(0) = 2^- The condition 
would be represented as 



P(0) 



C 



V2: 



■ sm- 



vr^ 



dy- 



y cosh I 



1 



(3.26) 



(cosh y + l)-y/(cosh y + cos b) 27rA 

The combination of ( |3.26| ), (|3.3| ) and (|3.9| ) provides the phase transition points from the 
lower gap to the two gap. By numerical calculations based on ( 3.26| ), (|3.3|) and (|3.9|) , we 
obtain the phase transition points plotted in Fig. ^(a) 



Let us see the phase transition points from the standpoint of the two gap phase. If we 
substitute a = into ( 3.16| ), it becomes the same as ( 3.26| ), which is the condition for the 
phase transition points. And if we set a = at ( |3.17 ), it becomes the same as (|3.9|). Let us 
check whether eigenvalue density in the two gap phase ( 3.1^ ) becomes the one in the lower 
gap phase ( |3.8| ) in a = limit. By using ( 3.26 ), we can replace /92,tg(A, 0, 6; a) of ( 3.1S| ) by 



c 



P2,t(;(A,0,6;a) = — ^ 



dy 



ycos ^ 



sm 



2 b 



,2 a 



(3.27) 



\/{l + cosh y) (cosh y + cos b) 
Then by summing up with Pi{„{C, 0, 6, c; a), ( 3. IS] ) at a = becomes (|3.8D . So we can see 



that solutions of two gap phase are smoothly connected to solutions of the lower gap phase 
at the phase transition point. 
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Figure 3: These are the plots of phase transition points in the regular fermion theory. Fig. 3(a) 
shows the plots of the phase transition points from the lower gap to the two gap phase, and Fig. 3(b) 
shows the ones between the upper gap and the two gap phase. 



Phase transition from the upper gap to the two gap In the upper gap phase, if 
the minimum of the eigenvalue density /9(vr) reaches to /j(vr) = the phase transition from 
the upper gap phase to the two gap phase will occur. The condition would be represented 
as 



C a 

^tCOS- 



dy 



ye 



(3.28) 



2 Jc "" v^(e-2j/ + 2e-ycosa + l)(l - e-y) 27rA 
The combination of ( p.28| ), ( |3.13D and ( |3.3| ) provides the phase transition points from the 
upper gap to the two gap. By numerical calculations based on ( |3.28|) , (|3.13| ) and ( |3.3D , we 



obtain the phase transition points plotted in Fig. p(b) . 

Let us see the phase transition points from the standpoint of the two gap phase. If 
we substitute 6 = vr into ( ^.16 ), it becomes (|3.28|) . And by using the relationship ( 3.2^ ), 
( 3.17 ) at 6 = vr is reduced to (|3.13| ). Let us check that the eigenvalue density in the two 
gap phase ( 3.18| ) at b = n is smoothly connected to the one in the upper gap phase ( 3.11| ). 
To confirm it, first we apply the following formula 



/i(a,7r, a) 



vr 



cos 



2 o 



sm^Wsm^l 



sm 



2 a 



COS ■ 



(3.29) 



and then by using ( 3.28 ) we can rewrite p2,tg as 



/02,is(A,a,7r;a) 



27rA 



c 



vr^ 



dy 



ye y sin ^ 



,2 a 



,2 a 



(3.30) 



(1 - e-y)\/l + 2 cos ae-y + e-^y ' 

Then by summing up with /9^'{ ((", a, vr, c, a), we can see that ( ^.18 ) at 6 = vr becomes 
( 3.11 ) . So we can see that solutions of two gap phase are smoothly connected to solutions 
of the upper gap phase. 

Quadruple phase transition point in the regular fermion theory There is a 
quadruple phase transition point 



a; 



rj 



0.596967, 



C^ = 2.864539029, (c^ = 0.644715) 



(3.31) 
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at which the no gap, the lower gap, the upper gap and the two gap phase coexist. Let us 
check whether (3.31) is the quadruple phase transition point. Please note that ( p.21[ ) is 
the condition for the phase transition from the no gap to the lower gap phase. And ( |3.2^ ) 
is the one between the no gap and the upper gap, ( p. 26 ) is the one between the lower gap 
and the two gap and ( 3.2g| ) is the one between the upper gap and the two gap. So if ( 3.31| ) 
satisfies these four equations simultaneously, it would be the quadruple phase transition 
point. If we substitute 6 = vr into ( |3.26| ) it becomes 



^ / dyy. 



1 



And if we substitute a = into ( 3.2^ ), it also becomes the same equation as (|3.32|) . We can 
also see that if both ( 3.21D and ( 3.24D are simultaneously satisfied, ( 3.32| ) is also satisfied. 



So if there is a point satisfying (3.32), the point becomes the quadruple phase transition 
point. By the study in |||], we already know that only the ( 3.31 ) simultaneously satisfies 
both condition ( p.21[ ) and ( 3.24 ). Hence ( p.31| ) satisfies ( |3.32| ). So ( 3.31 ) becomes the 
quadruple phase transition point in the regular fermion theory. Then due to the existence 
of the quadruple point, the phase structure of the regular fermion theory becomes as in 
Fig. 2(a) . We can also check that it becomes a = and 6 = vr at (|3.31 ), and then the 
eigenvalue density functions in each phase (|3.8D , (|3.1l|) ( ^3.18 ) and (6.33) of [j|] coincide. 



3.2 Critical boson theory 

In this subsection we will study the higher temperature phases, the lower gap, the upper 
gap and the two gap phase of the level k U{N) Chern-Simons matter theory coupled to 
massless critical bosons in the fundamental representation. The low temperature no gap 
phase is already studied in section 6.2 of |j|]. This is a Chern-Simons gauged version of 
the U{N) Wilson Fisher theory. The Lagrangian of the UV theory is that of massless 
minimally coupled fundamental bosons deformed by the interaction. 



5S 



d^xA<f>(f), 



(3.33) 



where A is a Lagrange multiplier field. ^ The effective potential in this theory is given by 
(3.12) of g and it is ^ 



ViU) 



N^C 3 ^ iv^c '•~ 



Qtt 



-a^ + 



2tt 



dy / da yp{a) (ln(l - £-2^+^") + ln(l - e-^^"*")) 



(3.35) 



=V^-''[p,N], 



''See subsection 4.3 of \^ for details, and in particular eq.(4.35) for the Lagrangian; note ||9j employs the 
symbol a for our field A. 



^Originally in 111], the effective potential is obtained as (3.9) of ||l|, 

ydy / p{a)da {[n{ 



N I ^ 2(cr^-A/?2)i\ iV 



1-e 



-y + ia 



) + ln(l - e-''-'-")) . (3.34) 



From this, constants A and a are obtained by requirement that they extremize the above (3.34) at constant 
p. By this, we yield A = a^T^, and the equation (p^3q). And by substituting these, the effective potential 
is simplified to be (3.35). 
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Figure 4: Phase diagram of the critical boson theory. Here (Ac,Cc) = (0.403033,4.24292) is the 
quadruple phase transition point where the four phases (no gap, lower gap, upper gap, two gap 
phases) coexist. 



where a provides the squared thermal mass as a T . And the a was determined by 

p{a) An 2 sinh(^^^) + In 2 sinh(^^i^)') = 0. (3.36) 

The general form of the free energy of the critical boson theory on S"^ x S^ is given as 

a-/3 



F^^ =y'-^[p, N] - N^V / da j dp p{a)p{p) log 
=V'-''[p,N] + F2[p,N]. 



2 sin ■ 



(3.37) 



To obtain the free energy in each phase, we only have to evaluate the eigenvalue density p 
and a in each phase and just substituting to ( |3.37D . Please note that a can be determined 
if the eigenvalue density is determined. So obtaining the eigenvalue density in each phase 
is equivalent to obtaining the free energy in each phase. 

For later use we will give the form of V'{z) here. We obtain V'(z) from ( p. 35 ), 



V'iz) = ^ I dy y 



-le 



le 



y-l 



e-y z - e~y 



(3.38) 



where z = e*". 



The phase structure of this theory is depicted in Fig. 4(a) . We will elaborate the phase 
structure in following subsubsections. 

3.2.1 Lovi^er gap phase 

Also in the lower gap phase of the critical boson theory, we employ the same procedure as 



the one in section 3.1.1 to obtain the eigenvalue density p. We employ the cut region and 



the cut function described in appendix AA. Because po = 0, U{z) is same as V'{z) 



U{z) = V'{z) = ^ r dyy 



-le 



y-l 



+ 



le 



z — e 



(3.39) 
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By substituting the above into ( 2.29| ), we immediately obtain H{u) as well as <5(n) as 



^u) = <^fg''{a,bX;u) 



'271 



dy y 



y^{e-y -e'^){e~y -e~ 

y roo 



-ib 



){u - e~y){u-^ - e-y) 



^2tt\,L '^^^u-e-y 



dy y- 



u 



e-y 



(3.40) 



From <I>+(n) — <1> {u) = Attp{u) at the cut u = e*" with —b < a < b, we obtain the 
eigenvalue density in the lower gap phase as 



p{a) 



C 



a 



sm 



V27r2 



sm 



dy 



y sinh | cos ^ 



-^cosh y — cos 6(cosh y — cos a) 



(3.41) 



And at the lower gap vr > |a| > 6, p{a) = 0. By substituing the eigenvalue density ( 3.41] ) 
into ( p. 361 ) and ( 3.37] ), we can obtain the free energy in the lower gap phase in the critical 
boson theory. 

The condition lim„_i.oo ^{u) = 1 requires following condition 



Mt^iC,<^,b)^-^ 



27r 



dx 



logx 



X 



log a; 



X Vx2 - 2x COS 6 + 1 



(3.42) 



By ( |3.42] ) and ( |3.36] ), we can obtain (6, a) as functions of (A, C) as (6, a) = (6(A, C), o"(A, ("))• 
( 3.41 ), ( 3.42] ) and ( 3.36] ) provide a set of complete solutions in the lower gap phase of the 
critical boson theory. 

3.2.2 Upper gap phase 

Next we will search for an upper gap solution in the critical boson theory. The cut region 
and the cut function h{u) in this case are defined in appendix A. 2] . Here the procedure to 
obtain the eigenvalue density function p is the same as the one in subsection p. 1.2 . We will 
take po{ct) as non-zero, and po{a) = ^^ at the upper gap region —a < a < a. 
Functions H{u) and ^{u) are obtained based on ( 2.29] ), 



^u)=^'^'g{a,a,C;u) + ^P,l{X,a;u 



Kq{(^, a, C; u) 



"27r 
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dy 
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ye-yh{u){l + e-y){l - u-'^] 
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^/{e-y - e^''){e-y - e' 

y roo 



'){u~^ - e-y){u - e-y) (3-43) 



dy 



ye 



271 \j„ "^ u — e y /„ " u ^ — e y 



dy 



ye~ 



Here we use the same definition of <l>(Jg(A,a; n) as ( 3.10] ). Prom (|3.43l) , by taking $"'"(n) — 
$~(tf) = 47r(/o(u) — po{u)) at the cut region u = e*" with vr > |a| > a, we can obtain the 
eigenvalue density function as 



p{a) 
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27rA V27r2 

c.b 
ug 



2 a ott 

sm sm — 

2 2 



dy 



y\ sin ^1 cosh | 



^cosh y — cos a (cosh y — cos a) (3.44) 



=P«g(C,A;<T,a;a). 
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In the region \a\ < \a\, p{a) = i^. By substituing the eigenvalue density ( |3.44|) into (|3.36| ) 
and ( p.37 ), we can obtain the free energy in the upper gap phase. To derive (|3.44| ), we 



have used the formula ( 3.12 ) 



From the condition, limu_5.oo ^{u) = 1 — /T dapo{a)^ we obtain the condition 

Here we have used the formula (|3.14| ). By (|3.45|) and ( |3.36|) , we can obtain {a, a) as 
functions of (A, C) as {a, a) = (a(A, C),cr(A, C))- ( ^-441) , ( |3.45| ) and ( ^.36 ) provide a set of 
complete solutions in the upper gap phase in the critical boson theory. 

3.2.3 Two gap phase 

We will search for a two gap solution in the critical boson theory. Our two cuts and the cut 
function are defined in appendix [A.3| . By the same procedure as the one in section 3.1.3, 
we can obtain H{u) as well as ^(u) as 

$t7 C, a, b, a- u) =^ dy ( \ _^ '-—- + ' \/l -— f 

'^ 2Tr J„ \{u — e y)(u ^ — e y) h[ey)[ey — u)[e y — u) 

(3.46) 

Here $4*0 (A, o, b; u) is already defined in (|3.15| ). From Imiu^oa ^{u) = 1 — f^^ da po{a), we 
obtain following two conditions, 



.T{a,b)=^y'^-'ia,b,a) 



and 



4vrA ' ' ' 27r' , _, 

y'-^a, b,a)= dy ^ =, 

J a y^(cosh y — cos a) (cosh y — cos b) 



^A(a,6)=l + ^g-V,M) 

\ (3.48) 



poo 

g--b{a,b,a)^ dyy 

J a 



ey + e-y 



Y^(cosh y — cos a) (cosh y — cos b) 



T{a,b) and A{a,b) are already defined in ( pB]) and ( ^ ). By using ( |3l8| ), ( |07D and 

( |3.36D , {a,b,a) are determined as functions of (A,C) as {a,b,a) = {a{XX),b{X,C),f^{^,C))- 

From ( ^3.461) and by taking $+(n) — <l>~(u) = 47r(/9(u) — /3o('u)) at the cuts \a\ < \a\ < \b\, 



we obtain the eigenvalue density function as 

P(") =Ptgioi) = p'i-^tg{C,a,b,c;a) + p2,tg{X,a,b;a), where 

Uc.bia,b;y) =^J {er'^y -2e-ycosa + l)(e-2y - 2e-s/cos6 + 1). 
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Here ( ^.49 ) shares the same definition of p^.tg-,^ with (|3.18D . p^^tg is the same functional 
form as the eigenvalue density in the GWW model (7.6) in |jl]. Pxiq can be regarded as 
an additional term depending on the detail of the theory. By substituing the eigenvalue 
density ( p. 49 ) into ( 3.36| ) and ( p. 37 ), we can obtain the free energy in the two gap phase 
of the critical boson theory. During the calculation of ( 3.49| ), we have used ( p. 47 ) also. 

The combination of ( p. 47 ), ( |3.48 ), ( 3.36| ) and ( |3.49| ) provides a complete set of solutions 
in the two gap phase. 



At large ^ limit, as we proved in appendix C.1.2 , the eigenvalue density approaches 
the universal distribution (p. 16 ) because a is kept as finite positive quantity at the limit. 
In the limit, a, 6 and the eigenvalue density behave as 



=7rA 



7rA + 



=8sin(7rA) exp 
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sin(7rA) 
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V 



cosh y — cos TT A 



(3.50) 
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■ cos 



a — a 
b — a 



At appendix C.2.2, we have demonstrated the above behavior ( p.50| ). As we can see in the 
appendix, the eigenvalue density is dominated by p2.tg at the large C limit. And the form 
of the eigenvalue density approaches cos~^ \/«i function same form as (7.11) in ^ which 
is the large C limit of the eigenvalue density of the GWW model. We can also see that 
the range of the domain of the cut e is smaller than the one in the GWW type matrix 
integration at the same value of (A, C)- 

3.2.4 Phase transition points 

Phase transition from the lower gap to the no gap Let us consider the behavior 
of lower gap solutions ( [3.42 ) and ( 3.41 ) at 6 = vr which would correspond to the phase 
transition points from the lower gap to the no gap. If we substitute 6 = vr into ( |3.42 ), it 
becomes 

T^ Jo 



1 



l + X 



n=l 



n^ 



(3.51) 



This is exactly same as the condition for the phase transition from the no gap to the lower 
gap phase which is obtained by substituting a = vr into eq. (6.24) of j|l| and by requiring 
p{tt) = 0. Based on this, if we substitute 6 = vr into ( 3.41 ), we obtain 
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y sinh | cos % 
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4(1 - e-y-^"^) 



(3.52) 
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n" 



In the last step we have used ( ^.51 ) and the Taylor expansion and an integration formula 
(3.23). The last line of ( p.52[ ) is exactly same as the eigenvalue distribution in the no gap 
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phase, described at (6.24) in Q. So, at the phase transition points, the lower gap solutions 
are smoothly connected to the the no gap phase solutions. 

Phase transition from the upper gap to the no gap Let us consider the behavior 
of upper gap solutions ( ^.44 ) and ( 3.45| ) at a = which would correspond to the phase 
transition points from the upper gap to the no gap phase. If we substitute a = into 



( 3.45 ), it becomes 
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vr ^^ 
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n^ 



(3.53) 



This is exactly same as the condition for the phase transition from the no gap to the upper 
gap phase discussed in |j|], which is obtained by substituting a = into (6.24) of ^ and 
by requiring p(0) = 2^. 

Based on this, if we substitute a = into (|3.44|) , we can see 



p{a) 



1 



c 



27rA 

1 
2^ 

1 

2^^ 



^/2^2 
7r2 



a 



sm 



c 



oo 

2^^ 

n=l 



dy y 



cos na 



sm 



1 



dy 



y\ sin yI cosh ■ 



"v/cosh y 
1 



cos a(cosh y — cos a) 



a=0 



2(1 - e-y 
1 + an 



4(1 



-y+ia 



) 4(1 - e-y-^'^) 



(3.54) 



n^ 



In the last step we have used ( ^.53 ) and the Taylor expansion and the formula ( 3.23| ). 
The last line of (|3.54D is exactly same as the eigenvalue distribution in the no gap phase, 
described at eq. (6.24) in [||. So, at the phase transition points, the upper gap solutions 
are smoothly connected to the no gap phase solutions. 

Phase transition from the lower gap to the t^vo gap In the lower gap phase at 
fixed A, if the maximum of the eigenvalue density p(0) reaches to p(0) = 2^^ ^^^ phase 
transition from the lower gap phase to the two gap phase will occur. The conditions would 
be represented in terms of ( |3.41| ) as 



P(0) 



C . h 
— 7: sm - 
27r2 2 



dy 



y^2(cosh y — cos b) sinh | 27rA 



(3.55) 



The combination of ( 3.55[ ), ( p. 36 ) and ( |3.42| ) provides the phase transition points from the 
lower gap to the two gap. By numerical calculations based on ( 3.55| ), ( |3.36 ) and ( 3.42 ), we 



obtain the phase transition points plotted in Fig. 5(a) 



Let us see the phase transition points from the standpoint of the two gap phase. If we 
substitute a = into (3.47), it becomes the same as ( |3.55| ), which is the condition of the 
phase transition points. And if we set a = at ( 3.481 ), it becomes the same as ( 3.42| ). Let 
us check that the eigenvalue density in the two gap phase ( 3.49| ) becomes the one in the 
lower gap phase (|3.41| ) in a = limit. By using (|3.55|) , we can replace /?2,tg(A,0, 6; q) of 



dMDby 



/52,tg(A,0,6;a) 



27r2 



dy 



y cos ^ 



sm^ I — sin^ ^ 



y^2(coshy — cos 6) sinh ■ 



(3.56) 
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Figure 5: These are the plots of phase transition point in the critical boson theory. Fig. 5(a)| shows 
the plots of the phase transition points from the lower gap to the two gap phase, and Fig. ^(b)| 
shows the ones between the upper gap and the two gap phase. 



Then by summing up with p^^-j ((",0, 6, c; a), ( ^.491) at a = becomes (|3.41 ). So, at the 
phase transition points, we can see that solutions of two gap phase are smoothly connected 
to the lower gap phase solutions. 

Phase transition from the upper gap to the two gap In the upper gap phase, if 
the minimum of the eigenvalue density p{tt) reaches to p{tt) = 0, the phase transition from 
the upper gap phase to two gap phase will occur. The condition would be represented as 



c 

—= COS - 

^/2^2 2 



dy 



y cosh I 



(3.57) 



2itX v^vr^ 2 7g. " ^cosh y — cos a (cosh y + 1) 
The combination of ( 3.571) , ( ^.45 ) and ( |3.36| ) provides the phase transition points from the 



upper gap to the two gap. By numerical calculations based on ( 3.57| ), ( |3.45 ) and ( 3.36| ), 



we obtain the phase transition points plotted in Fig. 5(b) . 

Let us see the phase transition points from the standpoint of the two gap phase. If we 
substitute b = tt into (|3.47| ) it becomes ( p. 57 ) which is the phase transition condition. And 
by using ( |3.57 ), we can see that ( |3.48D at 6 = vr is reduced to ( p. 451 ). Let us check that the 
eigenvalue density in the two gap phase ( |3.49| ) at 5 = tt is smoothly connected to the one 
in the upper gap phase ( p. 44 ). To confirm it, first we apply the formula (3.29). And then 
by using ( 3.57 ) we can rewrite p2,tg as 



P2,t<;(A,a,7r;a) 



27rA 
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ye ^|sin^ 



.2 a 
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(1 + e-y)Vl - 2cos ae-y + e'^v 



(3.58) 



Then by summing up with pfj. {(,a,Tr,c,a), we can see that ( |3.49|) at 6 = vr becomes 



( 3.44 ) . So we can see that solutions of two gap phase are smoothly connected to solutions 
of the upper gap phase. 

Quadruple phase transition point in the critical boson theory There is a quadru- 
ple phase transition point 
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0.403033, 



-C.6 



Cc = 4.24292, {ac = 0.644715), 



(3.59) 
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at which the four phases coexist. Let us check whether ( 3.59| ) is the quadruple phase 
transition point. Please note that (|3.51| ) is the condition for the phase transition from the 
no gap to the lower gap phase. And ( f3.53|) is the one between no gap and upper gap, ( 3.5^ ) 
is the one between lower gap and two gap and ( |3.57| ) is the one between upper gap and two 
gap. So if ( p. 59 ) satisfies these four equations simultaneously, it becomes the quadruple 
phase transition point. If we substitute 6 = vr into ( |3.55| ) it becomes 
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TT^ 



dy 



ye~ 
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1-e'^y 2-k\' 



(3.60) 



And if we substitute a = into ( 3.57| ), it also becomes the same equation as ( 3.60| ). We can 
also see that if both ( 3.51D and ( 3.53D are simultaneously satisfied, ( 3.60| ) is also satisfied. 
So if there is a point satisfying ( 3.60| ), the point becomes the quadruple phase transition 
point. By the study in |||], we already know that only the ( 3.59| ) simultaneously satisfies 
both condition ( |3.5lD and ( 3.53| ). Hence ( ^3.59 ) satisfies ( 3.60| ). So (|3.59 ) becomes the 
quadruple phase transition point in the critical boson theory. Then due to the existence 
of the quadruple point, the phase structure of the critical boson theory becomes as in 
Fig. 4(a) . We can also check that it becomes a = and 6 = vr at ( |3.59 ), and then the 
eigenvalue density functions in each phase ( |3.41| ), ( p.44[ ) ( |3.49| ) and (6.24) of |l|] coincide. 

3.3 Duality between the regular fermion theory and the critical boson theory 

In this subsection we will verify the level-rank duality between the level k U{N) critical 
boson theory and the level k U{k — N) regular fermion theory. We will show that the quan- 
tities at {Xc.b,Cc.b) in the critical boson theory have duality relationships to the quantities 



at (K.fXr.f) = (1 — K.b 



' i-A. 



-Cc.b) in the regular fermion theory as 
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(3.61) 



(3.62) 



br.f = vr — Oc.fe, a-r.f 

Here the p^.f is the eigenvalue density in the regular fermion theory and p^.b is the one in 
the critical boson theory. The relationship between {Xc.b-,Cc.b) and (K.fXr.f) is 
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c.b: 



Cr.f 
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i-x. 



-c 



c.bi 
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^ -X 
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'^r.f 



(3.63) 



If the above duality relationships are valid, the free energy at (Ac.6, Cc.b) in the critical 
boson theory agrees with the free energy at (A^./, Cr./) — \^ ~ ^c.b-, i_''x Cc.b ) in the regular 



fermion theory. So if we verify the above relationships ( p. 61 ) ~ ( 3.63| ), it would be the 
proof of the level-rank duality between the regular fermion theory and the critical boson 
theory. 

In this subsection, first, we will demonstrate the relationships ( 3.611) ~ ( 3.63 ) between 



the lower gap phase of the regular fermion theory and the upper gap phase of the critical 
boson theory. Next we will verify the relationships between the upper gap phase of the 
regular fermion theory and the lower gap phase of the critical boson theory. After that, we 
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will demonstrate the duality relationships between the two gap phase of the both theory. 
By these three steps, we complete the proof of the relationships ( 3.61 ) ~ ( |3.63D between 



the both theories. As the final part of this subsection, we will completely verify the duality 
by evaluating the free energy in the both theories. 

3.3.1 Between the low^er gap phase of regular fermion and upper gap phase of 
critical boson 

Let us define sets of quantities S^g,Sj^ as 

•^ug ={^c.b,Cc.b]ac.b,Cr,Pug{oi)), 



^c.b \2vrAc.b 



1 - K.b, z z—Cc.b] VT - Uc.b, cr, r — — Pug (vT + a, 



Please note that the eigenvalue density ( 3.44| ) and the two equations ( 3. 451) , ( |3.36| ) provide 



a complete set of solutions in the upper gap phase of the critical boson model, while 
the eigenvalue density (|3.8|) and the two equations (|3.9D , (|3.3D provide a complete set of 
solutions in the lower gap phase of the regular fermion theory. So we can verify the duality 
relationships by checking that S^^^ becomes a solution provided by the three equations 
in the regular fermion theory if S^'g is a solution provided by the three equations in the 
critical boson side. 

First we will confirm the relationship ( p. 611 ) between the eigenvalue densities ( 3.44| ) 



and ( |3.8| ) under the correspondence ( p.63| ) , ( p.62| ) . We can show the relationship by directly 
substituting the (|3^ ), (^) into (U) as ^ 
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^c.b f 1 _ c.b 
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(3.66) 

Pug iCc.b, K.b', O", «c.6; TT + a) 1 . 



Steps of the calculation to verify ( 3.66 ) are basically combinations of cos(2; + tt) = —cos a;, sin ■ 



cos ^ , as follows 
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1 — Ac.b \'2-k'^ V 2 2 j^ ^coshy + cos(7r — ac. b)(coshy + cos a) 

Ac.6 Cc.b L„2 7r-ac.6 „:„2 " f°° j,. y COS f cosh f 



sin TT / c^y 



1 — Ac.b \'2tt^ V 2 2 j^ ^coshy + cos(7r — ac.b){coshy + cos a) 

^pI'J I TI^T — Cc.6, 1 - Ac.b; cr, TT - Uc.b; a j . 

Basically, almost all steps of calculations to verify the duality relationships are just repeating the similar 
steps to the above. So from here we will not give a blow-by-blow description of similar steps to avoid to be 
lengthy 
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We can also see that S^^ becomes a solution of (|3.9|) if S'^^g is a solution of ( 3.45| ) by 
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(3.67) 



You can check this by a straightforward calculation similar to (|3.6 



We will also see that S^ becomes a solution of ( |3.3D if S^g is a solution of ( 3.3(]| ). 



This is already verified in section 3.6.2 of [|[ in the general case. We can check it as follows: 
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(3.68) 



Here we have used 



27r 



da (log(l + 2 cos ae"^ + e"^^) =0 (c > 0) 



(3.69) 



27r 



confirmed by the Taylor expansion and L dae*" = 



By the above discussion, we can conclude that S^' becomes a solution provided by 

and (|3.3D if 5^g is a solution provided by the equations 



the three equations (^;^^, ^tn^/ -"— - y ,^.^ ; ^- w^ic, 
( 3.44 ), ( 3.45 ), and ( p.36| ). Then we have finally verified the duality relationships between 
the lower gap phase of the regular fermion theory and the upper gap phase of the critical 
boson theory. 

3.3.2 Between upper gap phase of regular fermion and lower gap phase of crit- 
ical boson 

Let us define sets of quantities Sf;^, Sug as 



Sfg =iK.b,Cc.b;bc.b,cr,pfg (a)), 
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1 — Ac.b ■ ■ 1 — Xc.b V^TrAc.fe 

By the similar reason to the one in section ^.3.1 , we can verify the duality relationships by 
checking that Sug becomes a solution provided by the three equations ( 3.11| ) ( 3.13 ) and 



|3.3|) in the regular fermion theory if 5^^ is a solution provided by the three equations 



( 3.41D ,( 3.42 ) and ( 3.36| ) in the critical boson side. 

We will confirm the relationship ( 3.61| ) between the eigenvalue densities ( |3.41| ) and 
(3.11) under the correspondence ( 3.63| ),( p^6^ . By a similar calculation to ( p.65| ), we can 
show the relationship (|3.6l|) as 
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(3.71) 
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By checking 
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in a similar way to ( |3.67] ) , we can also see that Sug becomes a solution of ( |3.13D if Sf^'^ is 
a solution of ( 3.42| ). 

By repeating the analogous calculation to ( |3.68| ), we can see that Sug becomes a 



solution of (3.3) if Si^^ is a solution of (p. 3' 



So by the above discussion we can conclude that Sug becomes a solution of the three 
equations ( ^.H]) , ( |3.13|) and (^) if 5^^* is a solution of the equations (|3.41| ) , (|3.42D , and 



( 3.36 ). Then we have finally verified the duality relationships between the upper gap phase 
of the regular fermion theory and the lower gap phase of the critical boson theory. 



3.3.3 Between tw^o gap phases of regular fermion and critical boson theory 



Let us define sets of quantities S^^^Sl' as 
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By a similar reason to the one in section 3.3.1, we can verify the duality relationships by 



checking that Si' becomes a solution provided by the four equations ( p. 18 ), ( 3.16D ,( 3T^ ) 



and ( p. 3D in the regular fermion theory if 5f„ is a solution provided by the four equations 
(|3l9D,(|3^, (|3l^ ) and ( ^^361) in the critical boson side. 



Correspondence ( |3.61| ) betw^een the eigenvalue density functions First we will 
confirm the relationship ( |3.61| ) between the eigenvalue densities ( |3.49| ) and (3.18) under 
the correspondence ( p.63 ),( |3l62| ). We can directly check 
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Pi% ('^c.fe, Oc.b, bc.b, 0-; vr + a) (3.74) 



P2,tg (1 — Xc.b,ar.f,br.f', 



a) 



47r2(l - Xcb) A-/ h{uj) \u{uj-u)J ' 



(u + u) 



(3.76) 



^"We can see it from 

Ji (tt — 6, TT — a; TT + a) = — 2 



/•TI 

Jb 



d0 



'b (cos ^ — cos a) w (sin^ 
J^(7r — 6, TT — a; TT + a) =J-{a, b, a). 
These provide the complex hne integration along the lower gap 



|-sin2f)(sin^|-sin^|) 



(3.75) 
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P2,tg (Ac.6, VT - br.f, IT - a^./; TT + a) 



ih+{u) 



dbJ f {(jj + u) 



c.b Jl^J^ h{u}) \u{uj -u) 



(3.77) 



Here Vugs denotes the complex line integration along the upper gap running counterclock- 
wise from e"*"""-^ to e^""^-^ . And V^'^^ denotes the complex line integration along the lower gap 
running counterclockwise from e*^''-'' to e~'^^^f . By using the formula proved in appendix 



B.2.1 



du) u + u 



+ 



d(jj u + u 



-2m 



we can see that 



LuL ^(^) ^(^ ~ ^) Ar-f H^) ui^ - u) h+{u) ' 

^y^ Igs 

P2,tg (1 - Ac.b, TT - bc.b, vr - Ucb] «) 



(3.78) 



1 



A, 



c.b 



27r(l - A,.6) 1 - A, 



c.b 



P2,tg (Ac.fe, ac.b, bc.b] IT + a) 



(3.79) 



where we use a^,/ = tt — bc.b, K.f = vr — Ccfe. From ( [3. 74)) and (|3.79| ), we can confirm the 
relationship (|3.61| ) between the eigenvalue densities ( |3.49| ) and ( ^.181) , 



Pt/(a) 



A 



c.b 



1 — Ac.fe V 2-irX. 



1 



Vfe 



Ptgia + vr) 



(3.80) 



Relationship between ( p. 47] ) and ( |3.16 ) We will check that S^^ becomes a solution 



of ( p.l6| ) if Sfg'' is a solution of ( |3.47| ) . We can easily check 

y-^i-K - bc.b, TT - ttc.b, cr) = y^'^iacb, bc.b, cr). 
From ( 3.ig| ), by a similar calculation to ( 3.75 ), we can see 

duj „ , , , . . f du 



'^{ac.b,hc 



-4i 



J^c.b 

ugs 



h{uj) 



T(7r - 6c.fe,7r - Ucb) = 4i / 



Igs 



(3.81) 



(3.82) 



Here L'^gg denotes the complex line integration along the upper gap running counterclock- 



^ugs 



wise from e ^"'"■^ to e^"'"-'', and Lf^^ denotes the one along the lower gap running coun- 



^Igs 



terclockwise from e "■'' to e '=■''. From the formula proved in appendix B.2.2 
see 






Igs 



'^{O'c.b, bc.b) = T(7r - bc.b, vr - a, 



c.b) 



we can 



(3.83) 



Hence by taking into account Cc.b^c.b = (1 ~ ^c.b)Cr.f, we can conclude 

— - — T{ac.b,bc.b) = -^y'''^{ac.b,bc.b,<7) 
47rAc.b 27r 

=^ir~r^ ^ — ^^(^ ~ ^c.b, vr - CLc.b) = - — -, — i7-y-^(vr - bc.b, vr - Uc.b, cr). 

47r(l - Ac.bj 1 - Act, ^vr 



(3.84) 



This (3.84) shows that 5^„ becomes a solution of (3.16) if 5^''' is a solution of (p. 47) 
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Relationship between ( |3.48 ) and ( ^.171) We will verify that Si' becomes a solution 



of the ( |3.17| ) if 5f„ is a solution of ( |3.48| ). By a straightforward calculation, we can see 



Q''-f{'K - bc.b, vr - ttc.b, cr) = G^'^io-cb, bc.b, cr)- 
From ( 3.ig| ), by a similar calculation to ( 3.751) , we can see 



A(ac.fe, bc.b) = -4i / duj—--, A(7r - bc.b, tt - «£.&) = -4i / duj 
Jl?.!;, K^) Jl ■ 



Igs 



h{uj) 



By using the formula shown in appendix |B.2.5 

1 



1 



TTZ Tc.b 

UQS 



bJ 1 f ^ OJ 

n[u}) TTi J i^c.b n[u}) 

Igs 



(3.85) 



(3.86) 



(3.87) 



we can see 
1 



Cc.b, 



4ttX 



c.b 
1 



-A{ac.b, bc.b) = 1 + -j-G''' {ac.b, bc.b, o") 



47r(l - Xcb) 



Mt^ -bc.b, -n-- ac.b) 



Xr 



Cc.b, 



1 - Xcb 47r 



Q''--f (it- bc.b, TT- ac.b, (^) 



(3.88) 



This ( 3.88 ) shows that S^g becomes a solution of the ( 3.17 ) if Sig^ is a solution of (0.48 ). 

Correspondence between (3^) and ( [3.36 ) By repeating the analogous calculation to 
( |3.68 ), we can see that 5[„ becomes a solution of ( |3.3D if 5^^^ is a solution of ( [3.36 ). 

Summary of the proof So by the above discussion, we can conclude that Sl„ becomes a 
solution provided by the four equations ( 3.18| ), (|3.16| ), (|3.17|) and (|3.3|) in the regular fermion 
theory if S^^^ is a solution provided by the four equations (|3.49|) ,( ^^ ), ([3.48 ) and ( 3.36| ) 
at the critical boson side. Then we have finally verified the duality relationships at the two 
gap phase between the regular fermion theory and the critical boson theory. 

3.3.4 Between the phase transition points 

We have seen that the conditions for the phase transitions ( p. 211 ), ( p. 241 ), ( |3.55| ) etc, can be 
obtained as a certain limit of corresponding equations between which we have already seen 
the duality relationships. So we have already shown the duality relationship between the 
phase transition points, in fact. For example, ( 3.55 ) is obtained by ac.b -^ limit of the 
equation (|3.47] ) and ( 3.28 ) is obtained by br.f — )• vr limit of (|3.16|) . We have already shown 
that ( 3.16| ) maps to ( 3.47 ) under the duality, moreover we can see ac.b = T^—br.f = vr— vr = 0. 
Hence the phase transition points between the upper gap to two gap in the regular fermion 
maps to the one between the lower gap and the two gap phase in the critical boson theory. 
Now we will write the phase transition points between the A phase to B phase in the 
regular fermion theory as (A,B,RF), and the transition points in the critical boson theory 
as (A,B,CB). We can see following dual relationships represented by these symbols 



(No,Low,RF) ^ (No,Up,CB), (No,Up,RF) ^ (No,Low,CB), 
(Low,Two,RF) -^ (Up,Two,CB), (Up,Two,RF) «^ (Low,Two,CB). 



(3.89) 
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3.3.5 Free energy and completing the proof of the duahty 

Now we have proved that for any value of (Act, Cc.b)^ there are duahty relationships ( 3.61| )^( ^.63 ) 
between the critical boson theory and the regular fermion theory. 

Under these correspondences ( ^.611) ^( 13.631 ), we will show the free energy of the level 
k U{N) critical boson theory at {Xc.b,Cc.b) (see ( 3.37] )) agrees with the free energy of the 
level k U{k — N) regular fermion theory at (K.fXr.f) = (1 — ^c.b, jdx^Cc.b) (see ( |3.4D ), 



F^ = V'%'-\ N] + F2[p'-\ N] = V-f[p^-f, k-N] + F^ip'-'f, k-N]= F^^ 



r.f 



(3.90) 



where p"^'^ can be represented by p'^'^ through the relationship ( p. 61 ). 



As already shown in (3.37) of |l[, by a straightforward calculation similar to ( 3.6^ ), 
we can show V''-\p''-\ N] = V-f[p''-f, k - N]. 

Under ( |3.61| ) , F2[p^', k — N] is expanded in terms of p as 



F2[p''-^,k-N] =- N'^V 



da 



d(3 p^-\TT + a)p^-\TT + p)log 



2 sin 



a-f3 



+ 2N' 



N^ 



2TTXr 



V / da 



d^ P^-^t: + p)\og 



2 sin 



a- P 



1 



V da dp log 



,2vrAc.b 
By using the periodicity of p'^'^{a) = p'^-^{27r + a) and by using 

a — 13 



2 sin 



a — /3 



(3.91) 



V da loe 



2 sin ■ 



0, 



we can check 



So we can see 



F2[p'--f,k-N]=F2[p'-\N]. 



?N 



?k-N 



(3.92) 



(3.93) 



(3.94) 



-^ c.b ~ ^ r.f ' 

under the relationships ( |3.61| )^( |3.63 ). Now we have completed the proof of the level-rank 
duality between the level k U{N) critical boson theory and the level k U{k — N) regular 
fermion theory. 



4. Higher temperature phases of the Supersymmetric Chern-Simons mat- 
ter theory 

Next we will study the higher temperature phases of the M = 2 supersymmetric level k 
U{N) Chern-Simons matter theory with a single fundamental chiral multiplet. From here 
we will call the theory as "the SUSY CS matter theory" or "the SUSY theory" throughout 
this paper. The Lagrangian of the theory we study is presented in equation (2.1), (3.58) 
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of 11] (see also (6.1) and (6.20) of g). From (3.21) of Q, effective potential V{U) in this 
theory is obtained as 



ViU) 



6^|A| 



c — 6|A| / dap{a)Re / dyylogtanh 



y + ia 



(4.1) 



where the thermal mass (for both the boson as well as the fermion in the supermultiplet) is 
denoted by cT. The constant c above is determined by the requirement that it extremizes 
V{Uy, i.e. 

c + ia 



Re / daXp{a) log coth ■ 



(4.2) 



The effective potential can be interpreted as the sum of the regular fermion potential (| 
and the critical boson potential ( 3.351) . So the potential can be written as 



ViU) = Vr.fiU) + V,,t{U) 



(4.3) 



where 

Vr.fiU) 
Vc.b{U) 



N\ fc 



67r V A 



c'^ + 3 / dap{a) / dy y(ln(l + e"^"*") + ln(l + e 



—y+ia 



)) 



N'^C-s N'^C ''°° 



dir 



•-c^ + 



27r 



dy / da yp{a) (ln(l 



(4.4) 
,-y+ia^ + ln(l - e-2^"'°)) . (4.5) 



Hence it is convenient to utilize and combine the results in the previous section to obtain 
the eigenvalue densities and the free energy in the J\f = 2 supersymmetric case. The general 
form of the free energy of the supersymmetric CS matter theory on S'^ x S^ is given as 



FZ = V'^[p,N] + F,[p,N]. 



(4.6) 



As in the other CS matter theories, to obtain the free energy in each phase, we only have to 
evaluate the eigenvalue density p in each phase and substitute them into ( |4.6[ ) and (^) . So 
obtaining the eigenvalue densities in each phase is equivalent to obtaining the free energies 
in each phase. 

For later use we will give the form of V'{z) here. We obtain V'{z) = V^ Az) + V^^^{z) 
from (|4^) and (|^) as 






-le 



+ 



le 



2tt 



z + e-y z-^ + e-y 



-ie y ie y 

I r + - 



z — e 



-y 



(4.7) 



The phase structure of this theory is depicted in Fig. |6(a) , as obtained in following 
subsubsections. 
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Two gap 

Lower gap \ ^^'0 = ('^c. 4) 




(a) 

Figure 6: Phase diagram of the SUSY CS matter theory. Here (Ac,Cc) = (0.5,1.62509) is the 
quadruple phase transition point where the four phases (no gap, fower gap, upper gap, two gap 
phases) coexist. 



4.1 Lower gap phase 

Also in the lower gap phase of the SUSY CS matter theory, we employ the same procedure 
as the one in section |3.1.1| to obtain the eigenvalue density p. We use the cut region and the 
cut function defined in appendix |A.1| . Because po = 0, U{z) becomes the same as V'{z), 
we can immediately obtain H{u) as well as $(n) = h[u)H{u) by substituting the above 
( |47| ) into ( |2.29D . Since ( p.29| ) is hnear of V'{u), the $(n) is obtained as the sum of ( ^ ) 
and ( |3.4C1| ) with a = c. So <^(u) is obtained as 



^{u) 



cl>;/(c, 6, C; u) + <^%\c, b, C; u) = <^tg{c, b, C; u). (4.8) 

*" with —b < a < b, we obtain the 



Prom $+(w) — <I> {u) = 4:7rp{u) at the cut u = e 
eigenvalue density in the lower gap phase as 



/»(") = P?JiC, A; c, b; a) + p^^(C, A; c, 6; a) = pf^{C, A; c, b; a). (4.9) 

And p{a) = for vr > |a| > b. Also this is written by the sum of (3.8) and ( |3.41| ). By 



substituing the eigenvalue density ( [4.9| ) into ( [4. 21 ) and (4.6), we can obtain the free energy 
in the lower gap phase. 

The condition liniu^oo ^(u) = 1 requires following equation 



rr.f 



M,™(C, c, b) = M// (C, c, b) + Mf/(C, c, 6) = 1. 



(4.10) 



r-/ 



These M^'' and Mj^^^ are defined in 



I) and dg). 
By ( [4.10| ) and (|4.2|), we can obtain (6, c) as functions of (A, (") as (6, c) = (6(A, C), c(A, C))- 
( |4.9| ), ( 4.10 ) and ( [4. 2] ) provide a set of complete solutions in the lower gap phase of the 
SUSY theory. 

4.2 Upper gap phase 

To look for the upper gap solution in the current supersymmetric CS matter theory, we 



only have to follow the same procedure employed in sections 3.1.2 and 3.2.2. The upper 
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cut region as well as the cut function are described in appendix IA.4 The function H{u) 
as well as $(u) are obtained by using ( p. 291) , and the function ^{u) here is evaluated as 



$(u) = $2« (A, a; u) + $;/(c, a, (; u) + $^„^(c, a, C; u). 



~ug 



ug\ 



(4.11) 



Here ^ug{c,a,C',u) and ^u^g{X,a;u) are defined in ( |3.10D , ^'^g{c,a,C',u) is given in ( 3.43| ). 



From (4.11), by taking <l>^(n) — ^ (u) = 47r(p(ii) — po{u)) at the cut region u = e*" with 
71" > |a| > a, we obtain the eigenvalue density function as 



P(a) = PuqiC, A; c, a; a) + p'^^{C, A; c, a; a) 



1 



ug^ 



2tt\ 



pTqiC^^c.a^a), 



(4.12) 



.^•/ 



where p^t'/ and p'^g are defined in ( 3.11 ) and ( 3.441 ) respectively. In the region |a| < 



p{a) = 2^- By substituing the eigenvalue density ( 4. 121 ) into ( ^ ) and ( |4.6| ), we can obtain 
the free energy in the upper gap phase. 

From the condition, \\ra.u^oo ^{u) = 1 — /T ^0:^0(0^)5 we obtain the equation 



K/(C, c, a) + M^-'iC, c, a) = M-(C, c, a) 



"9 



«ff " 



(4.13) 



where Mug {C:C,a) and M^„''(C,c, a) are defined in ( p. 131 ) and ( |3.45 ) respectively. 

By ( 4.13| ) and ([4.2]), we obtain (a, c) as functions of (A, C) as {a,c) = (a(A, C), c(A,C))• 
(412), ( 4.13D and ( ^.2[ ) provide a complete set of solutions in the upper gap phase in the 
M = 2 supersymmetric CS matter theory. 



4.3 Twro gap phase 

We will search for a two gap phase solution. We only have to follow the same procedure as 
the one in section p.l.3| and p.2.3| . We will take the two cuts region and the cut function 



defined in appendix A. 3. By the same procedure as section 3.1.3, we obtain the functions 
H(u) and <l>(u) as 

^u) = </(C, a, b, c; u) + '^fg\C, a, b, c; u) + $^°(A, a, b; u). (4.14) 

Here ^Ig {(, a, b, c; u) and ^^^ (A, a, b; u) are defined in ( p.l5| ), and ^fg{(, a, b, c; u) is given 

( p6D . 



m 



From this ( |4.14D , by taking <l)+(n)-<l> (n) = 47r(p(n)-/)o(u.)) at the cuts |a| < \a\ < |6|, 
we obtain the eigenvalue density function as 



P(«) =Ptg{oi) = P2,tg{\ a, b; a) + p'^^lgiC, a, b, c; a) + pl%{C, «, b, c; a) 
=PiQ(C,A;a,&,c;a). 



(4.15) 



At 6 < \a\ < vr, p{a) = and at \a\ < a p{a) = 2^- P2,tg{X,o,,b;a) and p'A[C,,a,b,c;a) 



are defined in ( 3.1^ ) and pf^- {C,a,b,c;a) is defined in ( p.49| ). Similar to ( |3.18| ) and ( 3.49] ), 
( 4.15 ) can be interpreted as the sum of GWW type eigenvalue density p2^tg and the other 
additional terms depending on the detail of the theory. By substituing the eigenvalue 
density ( 1.15 ) into (4^) and ([1.6|), we can obtain the free energy in the two gap phase of 
the M = 2 supersymmetric CS matter theory. 
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From lim„ 



and 



^(u) = 1 — J_ da po{a), we obtain following two conditions, 
F"(a, b, c) =y'-f{a, b, c) + y'-\a, b, c), 



-l-A{a,b)=l + ^g'^{a,b,c) 
47rA 47r 



(4.16) 



g^"(a, b, c) =g''-\a, b, c) - g'-'fia, b, c). 



(4.17) 



Here T{a,b) and y^'^ are defined in ( |3.16| ), and y is in ( 3.47 ). And A(a,6) and g^'-' are 
defined in ( |3l7| ) and and g^^-^ is in (lOsI) . 

By using ([4.17| ), ( [4. 161) and (^), (a, 6, c) are determined as functions of {X,() as 
(a,6,c) = (a(A,C),&(A,C),c(A,C)). The combination of ( p^ ), (lilBD , (pj) and (|4l^ ) 
provide a complete set of solutions in the two gap phase of the A/" = 2 SUSY CS matter 
theory. 



At large ( limit, as we proved in appendix C.1.3 , the eigenvalue density approaches 
the universal distribution ( |2.16| ) because c is kept as finite positive quantity at the limit. 
In the limit, a,b and the eigenvalue density behave as 



a =ttX 



p{a) 



)sin(7rA) exp 
1 



6 = vrA + - 
sin(7rA) 



CA 



dy 2ycoshy 
cosh y — cos^ vrA 



+ 



(4.18) 



vr^A 



■ cos 



a — a 



At appendix C.2.3 , we have demonstrated the above behavior ( 4.18| ). As we can see in the 
appendix, the eigenvalue density function ( [4.15 ) is dominated by P2,tg- And the form of 
the eigenvalue density approaches cos~^ \/ai function same form as (7.11) in [^ which is 
the large Q limit of the eigenvalue density of the GWW model. We can also see that the 
range of the domain of the cut e is smaller than the one in the regular fermion theory as 
well as the one in the critical boson theory at the same values of (A, (")• 

4.4 Phase transition points 

Phase transition from the lower gap to the no gap Let us consider the behavior 



of lower gap equations (|4.10[) and (|4.9|) at 6 = tt which would correspond to the phase 
transition points from the lower gap to the no gap. If we substitute 6 = vr into ( 4.10|) , it 
becomes 



-2i. 



TT Jo 



dx 



logx 



2C 



vr 



n=0 



1 + c(2n + 1) 
(2n + l)2 



-{2n+l)c 



(4.19) 



This is exactly same as the condition for the phase transition from the no gap to the lower 
gap phase which is obtained by substituting a = vr into eq. (6.15) of |]l| and by requiring 
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/9(-7r) = 0. By using the condition ( [1.19D , if we substitute 6 = vr into ( [4. 91 ), we obtain 

PzT(C, A, c, vr; a) =i- + 1 £ cos(2n + l)a (^^^^f^e-(-^^)^) . (4.20) 

This is the same as the eigenvalue density function at the no gap phase of the supersym- 
metric CS matter theory, which is described at (6.15) of |||]. Also this can be written as 
the sum of ( |3.22| ) and ( |3.52| ). So the lower gap solutions are smoothly connected to the 
ones in the no gap phase at the phase transition points. 

Phase transition from the upper gap to the no gap Let us consider the behavior 
of upper gap solutions ( [4.13 ) and ( 4.12| ) at a = which would correspond to the phase 



transition points from the upper gap to the no gap phase. If we substitute a = into 



( |4.13D , it becomes 



oo 



This is equivalent to the condition for the phase transition from the no gap to the upper 
gap phase discussed in [j|], which is obtained by substituting a = into (6.15) of [|l[| and 
by requiring p(0) = g^- By using ( [4.21D , the eigenvalue density ( 4.12|) at a = becomes 



Pr,(C,A,c,0;a) =i- + -^ f cos(2n + l)a (^^^g^) e-^--^)^ (4.22) 

This is the same as the eigenvalue density function at the no gap phase of the SUSY CS 
matter theory, which is described by (6.15) of |Q]. So, at the phase transition points, the 
upper gap solutions are smoothly connected to the no gap phase solutions. 

Phase transition from the lower gap to the t^vo gap In the lower gap phase at 
fixed A, if the maximum of the eigenvalue density p(0) reaches to p{0) = g;^, the phase 
transition from the lower gap phase to the two gap phase will occur. The conditions would 
be represented in terms of ([4. 91 ) as 



1 \/^C f°° I y sin 2 cosh | y sin ^ sinh | 



A T^ Jc \ (cosh y + 1 ) y^ (cosh y + cos b) (coshy — l)y^(coshy — cos6) / 

(4.23) 
The combination of ( [4.23| ), ( [4.2[) and ( [4.10| ) provide the phase transition points from the 
lower gap to the two gap. By numerical calculations based on these ( [4.23 ), ( [4.2| ) and ( 4.10| ), 
we obtain the phase transition points plotted in Fig. 7(a)[ 



Let us see the phase transition points from the standpoint of the two gap phase. If 
we substitute a = into ( [4.16 ), it becomes ( [4.23| ). And if we set a = into (4.17), it 



becomes the same as (4.10|). Let us check whether the eigenvalue density in the two gap 



phase ( [4.15[ ) becomes the one in the lower gap phase ([4. 91) in a = limit. By using ( 4.23| ) 



and by similar calculations to ( 3.27 ) and ( [3.56| ), we can confirm ( [4.15 ) at a = becomes 



( |4.9| ) . So we can see that the solutions in the two gap phase are smoothly connected to the 
ones in the lower gap phase at the phase transition points. 
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Figure 7: These are the plots of phase transition points in the SUSY CS matter theory. Fig. 7(a 



shows the plots of the ones from the lower gap to the two gap, and Fig. [7(b)| is the ones from the 
upper gap to the two gap phase. 



Phase transition from the upper gap to the two gap In the upper gap phase at 
fixed A, if the minimum of the eigenvalue density p(7r) = />(— vr) reaches to /^(vr) = 0, the 
phase transition from the upper gap to the two gap phase will occur. The condition would 
be represented in terms of ( [4.12| ) as 



1 
A 



dy 



\f2y cos I sinh | 



+ 



\f2y cos I cosh | 



\/cosh y + cos a (cosh y — 1) \/cosh y — cos a (cosh y + 1) 



.(4.24) 



The combination of ([4.24|), ([4. 21) and ([4.13|) provide the phase transition points from the 



upper gap to the two gap. By numerical calculations based on (4.24), (|4.2|) and (4.13), we 



obtain the phase transition points plotted in Fig. [i"(b) . 

Let us see the phase transition points from the standpoint of the two gap phase. If 
we substitute 6 = vr into ( 4.161) and ( 4.17| ), they become ( 4.24|) and ( 4.13 ) respectively. By 



using these correspondence and by using the same procedure as the one in section 3.1.3 
and 3.2.3, we can also check that ( 4.15| ) at 6 = vr coincides with (4.12). So, at the phase 
transition points, the solutions in the two gap phase are smoothly connected to the upper 
gap solutions. 



Quadruple phase transition points in the SUSY CS matter theory There is a 
quadruple phase transition point 



ysusy 



0.5, 



/-susy 



1.62509, 



(c; 



-susy 
c 



0.542933), 



(4.25) 



at which the four phases coexist. Let us check whether ( 4.25| ) is the quadruple point. 
Please note that ( 4.231) is the equation for the phase transition from the lower gap to the 
two gap. ( 4.24| ) is the one from the upper gap to two gap, ([4.19|) is the one from the no gap 
to the lower gap, and ( 4.21| ) is the one from the no gap to the upper gap. So if the point 



(4.25) simultaneously satisfies the four equations, it becomes the quadruple point. We can 
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see that at a = 0, 6 = vr, ( 4.23 ) as well as ( [4.24 ) become 



2C 

2 



vr 



^T / dx 



logx 
1 — x"^ 



2ttX 



(4.26) 



We can also see that at A = 21 ( |4.19 ) and ( [4. 211 ) become equivalent to ( 4.26| ). By the 
study at [0, we already know that both ( 4.19] ) and ( f4.21| ) are simultaneously satisfied at 
the point ( 4.25| ). So also (4.26) is satisfied at ( |4.25p . Then we can see that (4.25) is the 
quadruple phase transition point in the SUSY CS matter theory. So due to the existence 
of the quadruple phase transition point, the phase structure becomes as in Fig. |6(a) . 



4.5 Self-duality of the SUSY CS matter theory 

We will study the self-duality of the supersymmetric CS matter theory under the Giveon- 
Kutasov type level-rank duality [||, |4|. The level k U{N) SUSY CS matter theory is 
dual to the same theory with level k U{k — N) gauge group. We will show that there are 
following relationships between the quantities in the level k U{N) theory at (A, C) and the 
quantities in the level k U{k — N) theory at (A,C) = (1 



A,T^C), 



A = 1-A, 



c 



A 



1-A 



c, 



-su / \ 

P (a) 



vr 



1- 

a, 



1 



A V27rA 



vr 



N _ k-N 



/"(vr + a) 



c. 



A 



(4.27) 

(4.28) 
(4.29) 



Quantities in the U{k — N) theory are denoted with over-line while the quantities in the 
U{N) theory are described without over-line. First we will see the duality between the 
quantities in the lower gap phase of the U{N) theory and the ones in the upper gap phase 
of the U{k — N) theory. Next we will demonstrate the duality in the two gap phase 
between the U{N) and the U{k — N) theories. At the end of this subsection, we will see 
the agreement of the free energy under the duality. 

4.5.1 Between the lower gap phase of the U{N) theory and the upper gap phase 

of the U{k - N) theory 

Let us define sets of quantities 



osu 

^19 



csu 



=(A,C;6, 



c,0 



1-A, 



1-A 



C;vr 



b,c, 



1-A V 27rA 



pf^iTT + a] 



(4.30) 



Please note that the the eigenvalue density ( ^^ ) and the two equations ( |4.1(]| ), (4.2) provide 
a complete set of solutions in the lower gap phase, while the eigenvalue density ( 4.12| ) 
and ( 4.131) , (|4.2| ) provide the one in the upper gap phase. So we can verify the duality 
relationships by checking that 5^" becomes a solution provided by the corresponding three 
equations in the upper gap phase if Sj^J' is a solution provided by the other three equations 
in the lower gap. 
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First we will confirm the relationship ( 4.28| ) between the eigenvalue densities ( [4.9D and 
(4.12) under the correspondence ( [4.27 ),( 4.29| ). Here by utilizing relationship ( |3.71| ) and 



( 3.66 ), we can check 



Puq(C,A;c,7r-6;Q) 



A 



1 - A V 27rA 



1 



pf«"(C,A;c,&;7r + a) . 



(4.31) 



This shows the relationship ( 4.28| ) between the eigenvalue densities (4.9) and ( 4.12| ). 

We can also see that S^ becomes a solution of ([4.13|) if 5f" is a solution of (|4.10| ) by 

showing 

1 — A -,„„ / A 



M4"(C,c,fe) 



A 



-M!, 



"n 1 - A 



C,c,7r 



To show this we have used the combination of ( 3.67| ) and ( |3.72| ) 



We will also see that the SfJt becomes a solution of the 

ay 



(4.32) 



if Sf^ is a solution of it. 

'9 



By using (4.31), we can directly check as 



Re / daA/9f"(a)logcoth 



c + ia 



-^ c 



<^ c 



Here we have used 



da log coth 



-2 
-2 

X + ia 



Re / da{l - A)p*"(7r + q) log coth 



c + ia 



Re / da{l - X)p'^„ia) log coth 



c + la 



0, coth( 



X + iTT , 



X 

tanh(-) 



coth(-) 



(4.33) 



(4.34) 



( 4.33 ) shows that S^^ becomes a solution of ( |4.2| ) if 5f" is a solution of it. 

By the above discussion, we can conclude that 5*" becomes a upper gap solution 
provided by the three equations ( [4. 121) , (|4.13D and (p^ at U{k - N) theory if Sig is a 
lower gap solution provided by the equations (O), ( |4.10D , and (|4.2|) at U{N) theory. Then 
we have finally verified the self-duality relationships between the lower gap phase of the 
level k U{N) SUSY CS matter theory and the upper gap phase of the level k U{k — N) 
SUSY CS matter theory. 



4.5.2 Self-duality in the two gap phase 

Let us define sets of quantities 



5-=(A,C;6,c,p-(a)), 
{X,C,d,b,c,pf^{a)) 
A 



osu 



1-A,. 



C; TT — 5, vr — a, c, 



A 



1 



PtT (^ + «; 



(4.35) 



1-A" ' '1-A V2vrA 

Here Ptq{a) is originally defined as the eigenvalue density function of the dual variable 

pf^iC, A; o, b, c; a) = pf^iC, A; a, b,c;a + vr). (4.36) 
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Please note that the the eigenvalue density ( [4.15 ) and the three equations ( 4.16 ), ( 4.17| ) 
and (|4.2D provide a complete set of solutions in the two gap phase. So we can verify the 
duality relationships by checking following: If Sf^ is a two gap phase solution of the U{N) 
theory provided by the four equations, also Sfg becomes a two gap solution of the U{k — N) 
theory provided by the same equations. 



Correspondence ([4.28|) between the eigenvalue density functions First we will 



confirm the relationship ( 4.28| ) between the eigenvalue densities pfi and p|". By using 
( 3.74 ) and ( 3.7S| ), we can see 



pf"(C,A;a,6,c;a)=/>?"( 



A 



^ ^ 1 - A 
1 



C, 1 — A; vr — 6, vr — a, c; a) 
A 



(4.37) 



27r(l - A) 1 - A 



Ptg(C,A;a,6,c;7r + Q). 



So there is a relationship ([4.28| ) between pf" and pf" 



About (4.16) Let us check Sfg becomes a solution of ( [4.16 ) if Sfg is a solution of it. By 
using XQ = {1 — A)C and ( p. 83 ), ( 3.81[ ) we can see 



47rA 



T(a,6) 



2-K 



y'''(a,b,c) 



-T(7r — b,-K — a) 



A C 



y'''{Tr-b,Tr-a,c). 



(4.38) 



47r(l-A) '■■ "'■' ' 1-X2tt^ 
This shows that Sfg becomes a solution of ([1.16[) if 5f" is a solution of it. 



About ([4.17[) Let us check 5f^ becomes a solution of ( [4.17 ) if S^g is a solution of it. By 
using ( 3.85| ), ( [3.86[ ) and ( [3.87 ), we can see 



47rA 



A(a,6) = l + -^g^"(a,6,c) 
47r 



1 



47r(l - A) 



A(7r - 6,7r- a) = 1 + 



A C 
1-A47r 



g^"(7r-6,7r-a,c). 



(4.39) 



Please note that 



g^"(7r-6,7r-a,c) = -g""(a,6,c). 



(4.40) 



This ( 4.39 ) shows that Sfg becomes a solution of the ([4.17] ) if Sf^ is a solution of it 



About ([4.2]) By repeating the analogous calculation to ( 4.33| ), we can see that Sfg be- 
comes a solution of the (4.2) if 5f" is a solution of it. 



Summary So by the above discussion we can conclude that 5*" becomes a solution 
provided by the four equations (p^),(p^,(p7[) and (|3) in the level k U{k-N) SUSY 
CS matter theory if S^ is a solution of the same equations in the the level k U{N) SUSY 
CS matter theory. Then we have finally verified the self-duality relationships at the two 
gap phase between the level k U{N) SUSY CS matter theory and the same theory with 
level k U{k — N) gauge group. 
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4.5.3 Between the phase transition points 

We have seen that the conditions for the phase transitions ( 4.19 ), ( 4.21 ), ( 4.23| ) and (|4.24| ) 
can be obtained as a certain limit of corresponding equations between which we have 



aheady verified the duality relationships. So, by the same logic as in section 3.3.4 , we have 
already shown the duality relationships between the phase transition points, in fact. By 
using the analogous symbols to the ones used in (|3.8g| ), we can see the following duality 
relationships 

(No,Low,SUSY) ^ (No,Up,SUSY), (Low,Two,SUSY) ^ (Up,Two,SUSY). (4.41) 

4.5.4 Free energy and completing the proof of the duality 

Now we have proved that there are self-duality relationships ( 4.27| )^( 4.2g| ) in the SUSY 



CS matter theory. Under these correspondences ( 4.27 )^( |4.29D , we will show that the free 
energy of the level k U{N) SUSY CS matter theory at (A, C) agrees with the free energy 

A, -ArC), namely, 



of the same theory with level k U{k — N) gauge group at (1 — xv, ^^^ 



F: 



k-N 



y""[p"", k-N]+ F2[^"", k-N] = y""[p"", N] + F2[p"", N] = F, 



N 
sw 



(4.42) 



Here we have used ([4.61) . As already shown in (3.30) of [|[, by a straightforward calculation 
similar to ( |4.33| ), we can show y"[p^",fc - iV] = y"[p'*",Ar]. We can also check that 
F2[p, k — N] = F2[p, N] by the almost same calculation as ( p. 911 ) and ( |3.92| ). 
So we can see 



F. 



N 



F: 



k-N 



(4.43) 



under the relationships ( [4.27| )^ p.29p . Now we have completed the proof of the level-rank 
self-duality between the level k U{N) SUSY CS matter theory and the level k U{k — N) 
SUSY CS matter theory. 



5. Analytic proof of the duality in the two gap phase of the GWW type 
matrix integral 

In [j|], phase structure of the toy large A^ GWW type matrix integral has been already 
investigated. Although several numerical evidence have been provided, the analytic proof 
of the level-rank duality in the two gap phase of the GWW model has not been completed. 
Then we will complete the analytic proof of the duality by using the three formulae ( 3.7^ ), 
dPI) and (|3^ ). 

This proof would be also useful for analysis of general CS matter theories on S'^ x S^. 
As we can see at ( 3.20[ ), ( p.50|) , ( f4.18| ), and appendix C^, the eigenvalue densities converge 
to the one in the GWW model at large ( limit. Moreover, every eigenvalue density of CS 
matter theories at the two gap phase includes the term p2^tg which is the same functional 
form as the eigenvalue density in the GWW model. So the duality relationships shown 
in this section play one of the fundamental role to establish the level-rank duality in CS 
matter theories on 5^ x S^. 

Let us focus on the proof of the level-rank duality in the two gap phase of the GWW 
model. As described in section C.4 of [P, we can regard (A, () as functions of the edge of 
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the cuts (a, b). To complete the proof of the self-duahty of the GWW type model in terms 
of this, we should show followings: 

A(a, b)({a, b) = A(7r — a, vr — 6)C(vr — a, tt — 6), (5-1) 

A(7r - 6, vr - a) = 1 - A(a, b), (5.2) 

Ptg f (1- A),^— ^,aj = J— J (^-^ -ptg{XX,a + -^)] ■ (5.3) 

If we can verify the above three equations, it becomes the proof of the level-rank duality 
in the two gap phase of the GWW type matrix integral. 



5.1 Proof of 

From (C.17) of |l| and by identifying w = e**^, we can see 



i-JT It hiuj) 

From this, by a direct substitution, we can also see 



C{a,b)X{a,b) = - f -^. (5.4) 



1 r ri 

C{7r-b,TT-a)X{7r-b,7r-a) = -— ——. (5.5) 



So by using ( [3.83 ), we can immediately verify ( ^.ij) . 



5.2 Proof of (^) 

From (C.17) and by the identification u = e , we can see that A(a, b) can be written as 

X(a, b) = — / -— — - I uj (cos a + cos 6) ) . (5.6) 

By the same way, we can also see 

A(7r — b,TT — a) = — / -— — - ( u) (cos a + cos 6) I . (5-7) 



So by using both ( 3.83| ) and ( |3.87| ), we can show that 

A(7r-6,7r-a) + A(a,6) = 1. (5.8) 

This is the proof of (|5.2|). 



5.3 Proof of (^) 

The eigenvalue density (7.6) of |j^] can be described by the complex line integral as 

p{X{a, b), C(a, b),a) = '-^ [ du; -]- f ^- + -) , (5.9) 
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here we identify uj = e . In the same way, by using (|5.2| ) and (5J.), we can also check that 

ih'^{u) f doj 



p {X{tt — b,TT — a), C(vr — 5, vr — a), vr + q) 



2 1 

+ - 



47r2(l - A) Jl^^^ h{uj) \uj-u u 
So by using ( p. 78 ), (5.2) and ( ^.1[ ), we can check 



p{X, C, a) -\ ^ — p ( 1 - A, :; ^C; T^ + a 



A 



1-A 



27rA 



(5.10) 



(5.11) 



This is the proof of (5.3). 

Then we have verified (5.1),( |5l^ ) and ( |5.3D . Hence we have completed the analytic 
proof of the level-rank duality in the two gap phase of the GWW type model. 



6. Summary and discussions 

In this paper, we have elaborated the higher temperature phases of the Chern-Simons 
matter theories with fundamental representation on S'^ x S^. Here we studied the lower 
gap, the upper gap and the two gap phases of the regular fermion theory, the critical boson 
theory, and the M = 2 super symmetric Chern-Simons matter theory. We have obtained the 
eigenvalue density function and the free energy in each phase, and we have also obtained 
the phase transition points and the phase diagrams. Based on these obtained quantities, we 
have completely analytically verified the level-rank duality (Giveon-Kutasov type duality) 
between the regular fermion theory and the critical boson theory, and the self-duality of 
the SUSY CS matter theory. We have also supplied the analytic proof of the duality in the 
two gap phase of the GWW type matrix integration which is suggested in Q. 

We can see that the proof of the duality in the two gap phase of the GWW model as 



well as the three formulae ( p. 78]) , ( 3. 83 ) and ( p.87|) proved in appendix B^ are based on 
the cut structure of the two gap phase, which would be shared by all CS matter theories 
on S"^ X S^. Please note that p2,tg in the eigenvalue density functions at the two gap phase 
is the same functional form as the eigenvalue density of the GWW model, and it is shared 
by the all eigenvalue densities in the two gap phase of the all CS matter theories dealt 
with in this paper. And the three formulae play the role to establish the level-rank duality 
with respect to the P2,tg- Hence the duality structure of the GWW model and the three 
formulae would play the most fundamental role in the level-rank duality of general CS 
matter theories on S^ x S^. 

It is also interesting to compare the three CS matter theories in this paper. As we have 
seen at ( 3.20 ), ( 3.50D , ( [1.18 ), and appendix p.2| , in the supersymmetric theory, the range of 
the domain of the cut e is smaller than the one in the regular fermion theory as well as the 
one in the critical boson theory at the same value of (A, C). And we can see that the phase 
transitions at the SUSY theory occur in lower temperature than the ones in the other two 
CS matter theories. (See Fig. 8(a) , g(b)| .) From this we can guess that there should be 
stronger attractive force between the eigenvalues in the SUSY theory than the one in the 
other two theories. This would be because the SUSY theory has more degree of freedom 
of matter fields than the other two theories. And then the effective potential in the SUSY 
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VjT- 



V^T- 





(a) 



(b) 



Figure 8: Fig. 8(a) is the phase transition points of the SUSY theory v.s the ones in the regular 
fermion theory. And Fig. p(b)| is the phase transition points of the SUSY theory v.s the ones in the 
critical boson theory. Here bold blue line in the both figure is the plots in the SUSY theory. The 
red dashed line in Fig. j(a)| is the points in the regular fermion theory and the black dashed line in 
in Fig. S(b) is the points in the critical boson theory. From these we can see the phase transition 



points in the SUSY theory locate in lower temperature than the other two theories. 



theory obtained after integrating out the matter fields, which causes the attractive force 
between the eigenvalues, are written by the sum of the ones of the regular fermion and the 
critical boson theory. Hence the attractive force would be enlarged by the summation. 

At the very high temperature limit, the eigenvalue densities of these three theories 
and the one of the GWW model converge to cos~^ ^/a, and they finally become the same 
universal distribution (|2.16| ) if we further increase the temperature. We have given the 
analysis in (|3.20|) , (|3.5C1| ) , ([4.18| ) and appendix ^. So at the very high temperature limit the 
results coincide with the ones in Q, we can regard their results as the very high temperature 
limit of our results. Our result can be regarded as a lower temperature generalization of 

i- 

It is easy and straightforward to study the critical fermion and the regular boson theory 
on S'^ X S^ in the same way as this paper. 

One of the interesting application of this paper would be to analyze the phase struc- 
ture of the dual gravity description governed by parity violating Vasiliev's higher spin 
equations |3|, §, ^ at finite values of the 't Hooft coupling A = -^ |^ 0]. The thermal 
system studied in this paper is dual to the Euclidean Vasiliev system in global AdS space 
compactified on a circle on circumference 2itR = y. The field theory analysis performed 
in this paper implies that this Vasiliev system admits a classical limit in large A^ limit 
with V2T^ = C^ aiid C held fixed. The saddle points obtained in this paper should map to 
classical solutions of this bulk description; however the equations of motion governing this 
classical description are not necessarily Vasiliev's equations. This is because quantum cor- 
rections to Vasiliev's equations proportional ^^^y— are not suppressed compared to classical 
effects in the combined high temperature and large N limit studied in this paper; such 
terms could modify Vasiliev's classical equations. It would be fascinating to determine the 
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(a) 



(b) 



(c) 



Figure 9: Fig. 9(a)| , |9(b)| , |9(c) show the cut region in the lower gap phase, in upper gap phase, in 
the two gap phase, respectively. 



effective 3 dimensional bulk equations dual to the large N limit of this paper, and to study 
the classical solutions dual to the saddle points described in this paper. 
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A. Cut function h{u) 

We will list the cut region and the cut function in this appendix. The cut region in each 



phase is depicted in Fig. 9(a) , 9(b) and 9(c) 
A.l Lower gap case 



Jb 



In the lower gap phase, there is a single cut from Ai = e , to Bi = e counterclockwise. 



As depicted in figure |9(a)| , the center of the cut is located at the point 1. 

Now let us define the cut function h{u) carefully. Let u = \u\e with 9 G (— -7r,7r], then 
we will define the h{u) as 



h{u) = ({u-e"'){u 



1 

-ib\\ 2 



(A.l) 



The single cut of this function is taken to lie in an arc on the unit circle extending from 
Ai to Bi. The function h is discontinuous on this arc with 



h^iz) 



iB_ I 96 r. 9 

±2e 2 A/ sm^ sm^ - 

2 2 



(A.2) 



46 



where z = e . Away from this cut (i.e. on the gap) on the unit circle 



ie 



/i(z) = 2ie'5'^/sin2|-sin2| e>0, 



h{z) = -2ie^ y'^sin^ | - sin^ | e<0. (A.3) 

Recah that sin^ ^ — sin^ 2 everywhere on the gap. On the real axis h{x) is everywhere real 
and is positive for x > 1, but is negative for x < 1. In particular h{0) = —1. In more detail 
along the real axis 



h{x) = \/ {x — cos b) + sin 6 > x > 1, 



h{x) = -J{x-cosb) +sm^b<0 x < 1. (A.4) 

A. 2 Upper gap case 

Let us consider the case with one upper gap without no lower gap. In this case there is 
a single cut extending from Ai = e*" to Bi = e~*° counterclockwise along the unit circle. 



The cut centers on the point -1 in the complex plane, see Fig. 9(b) . We can also regard 
this region as the complement of the cut in the lower gap phase. This interpretation would 
have important role to show the duality later. 

In this case, the function h{z) = y^(e*" — 2;)(e~*" — z) has a branch cut running coun- 
terclockwise from e*" to e~*°. The function h{u) obeys followings: 



-/i-(e^") = /i+(e^")= + 2ie*t,/sin2--sin2- if tt > a > a, 



2ie 2 . /sin — — sin — it — vr < a < —a, 



/i(e*") =2e*2 ^/sin^- -sin^ — if a>a>-a, 



(A.5) 



/i(x) = + y/(3; — cosa)2 + sin^ a for x > 1, 



2; — cosa)^ + sin^ a for x < — 1, 



h{z — )• 00) — T-z, 
h{0) =1. 

A.3 Two gap case 

In the two gap phase, our two cuts extend along unit circle from Ai = e*" to Bi = e , 
and from A2 = e~ to B2 = e~*" counterclockwise respectively, see the Fig. |9(c)| . The cut 
function h{u) in this two gap case will be 

h{u) = ((n - e")(u - e-'"'){u - e'^){u - e"*^)) ' , (A.6) 

where this function has branch cuts along the arcs enumerated above. We will summarize 
some of properties of the analytic function h{u). Along the unit circle outside the cuts, h 
is analytic as 
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/i(e*'^)=4e*%/(sm2--sin2-)(sin2--sm2-), -a < 6 < a, (A.7) 



/i(e*^) = -4e^Y(sin2--sm2-)(sm2--sin2-), -tt < 6 < -b, b < 9 < tt. (A. 
Along the cuts, h becomes discontinuous with 



h^{e'') = ± 4ie^%/(sin2 ^ - sin^ |)(sin2 ^ - sin^ ^), a < 9 < b, (A.9) 



/i±(e^^)=^4ie'V(sin2--sin2-)(sin2--sin2-), -b < 9 < -a, (A.IO) 

where /i = /i+ on the outside of the circle and h = h^ on the inside. Moreover h is real 
and positive along the real axis; explicitly for real x. 



h{x) =\/{l — 2xcosa + x'^){l — 2xcos6 + x^), — oo < x < oo. (A. 11) 

Note in particular that h{0) = 1. At large u we have 

/i(m) ~ +M — M(cosa + cos6), u ~ cxd. (A. 12) 

B. Important formula 

B.l Important formula in the upper gap phase 

In this subsection, we will derive several formulae in the upper gap phase with cut function 
h{z) s.t. h{z) = \J{e^^ — z)(er'^°- — z). This has a branch cut running counterclockwise 
from e*" to e~", and has the property ( |A.5| ). 



B.1.1 Proof of (|1|) 



Here we will prove the formula ( 3.14 ). 

hie except the cuts, from the Cauchy's theore 

(B.l) 



Since w^ is everywhere holomorphic except the cuts, from the Cauchy's theorem 



Please note from h{uj) — )■ w at |a;| — >■ oo, 



then 



r du 


r duj 

Jccuts K^) 


jj\ — >■ oo. 




/ doj 

J\lo\=oo K^) 


= 2m, 


d. _r 


d. _r 



(B.2) 



dw 

2m 



Cents H^) 71^1=1+161 H^) 7|w|=l-|e| hC^) 

duj f du! 

'\oo\ = l+\^\ K^) J\u\ = l-\e\ H^) 

-2 [ ^. 

J Lugs H^) 
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From first to second line we have used the fact that j^j^ is holomorphic inside the circle 

cond to third, we have u; 
Then we have finally derived 



wl = 1 — lei. And from the second to third, we have used ,.\ •, = — ,\ ■, in the cut region. 

Ill ' h+{ujj h (ui) *= 



J Lugs H^) 



B.1.2 Proof of (|3l^) 

Here we will prove the formula ( 3.12| ). 



Since , , w _ s is everywhere holomorphic except the cut, from the Cauchy's theorem 



duj f duj 



|<^|=oo h{uj){uj - u) Jc^^^^ h{uj){uj - u) ' 



(B.5) 



From TTTT 7 -^ ^ ^5 at |a;| — )• cxd. 



/|a;l=oo h{uj){uj-u) 



T\uj\=oo h{uj){uj-u) 

SO 

f doj f dw f dw 

J Cuts h{^){^ - U) J\oo\ = l+\e\ h{uj){uj - U) J\uj\ = l-\e\ h{uj){uj - u) 

duj f duj 

+ 



\u,\ = l+\e\ HUJ){UJ - U) J\uj\ = l-\e\ h{uj){uj - u) 

r duj r du / 1 1 

J Lugs K^){^ - u) Jl,^^^ h+{u) \{u-u)+ [u-uY 

(B.7) 

Here Lares denotes the complex line integration counterclockwise along the cuts. Please 
note on the cut region, ^^ 

^ --iTTiSioj-u), (B.8) 



{ijJ — ti)+ {ijJ — u) 

so it becomes, 

duj f I 1 \ _ 27rf 

Then we have obtained 

duj Txi 

,„^^ h{u:){u-u) = ~h^y ^^-^"^ 

B.2 Important formulae for the proof of level-rank duality in the tvvro gap phase 

In this subsection, we will derive several formulae in the two gap phase with cut function 

h{u) s.t. 

h{u) = [{u - e'''){u - e-'^Yu - (^^){u - e"^^)) ' , (B.ll) 



^^We should note from the standpoint of the counterclockwise complex integration, the sign of the delta 
function is flipped due to the opposite direction of the integration measure direction. 
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where this function has branch cuts extending along unit circle from Ai 



e'" to B^ 



e and from A2 = e~ to B2 = e~*" running counterclockwise. And this function has 
properties (A. 7) ~ ( A.12| ). Please note that the formulae derived here are important for 
the proof of the duality in the two gap phase. 



B.2.1 Proof of ( P78| ) 
We will show (|]7|), 



ih+{u) 



dcj Ficj) + 



J^un 



doj F{u}) 



2-k\ 



where we define F{uj) as 



Fiu) 



1 



2 1 

+ 



(B.12) 



(B.13) 



h{uj) \ijj — u u ^ 

Here throughout this subsection Ligg is the arc running counterclockwise along the lower 
gap on the unit circle, i.e, arc e*^ s.t. 6 < 16*1 < vr. And L^gs is the arc running counter- 
clockwise along the upper gap on the unit circle, i.e, arc e*^ s.t. —a < 9 < a. 

Since F{uj) is holomorphic everywhere except the cut region, due to the Cauchy's 
theorem, it becomes 

= (/ duj F(uj) - i doj F(oj). (B.14) 



doj F{uj) - <h doj F{uj). 

|a;|=oo J Cents 



From F{u!) — )• ^ at \uj\ — )• 00, we can see 



(p du F{uj) = j) duj F{uj). 

J\u]\=QO J Cents 



(B.15) 



From this, we can rewrite 



= ^ did F{uj) 

' ^cuts 



dhj 



\u\ = l+\e\ 



du} 



F{lj) - <P duj F{o 

JH = l-\e\ 

F{u;) + I duo F{u 

J\L0\ = l-\e\ 



M=i+k| J\^\=^-\A 

du {F+{u:)+F-{u:)). 



(B.16) 



From the first line to second line we have used the fact that F{(jj) is holomorphic inside 
the circle Iwl = 1 — lei. 



duj F{uj) = 



'|a;|=l-|e| 

Please note at the cuts region, 

1 



® doj F{ijj). 



(B.17) 



F+H + F-(, 



UJ 



h^{uj) \{uj — u)^ {oj — u) 



h+iu] 



{iTri6iu - u)) , (B.18) 
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and F^{uj) = F (u) = F{uj) in the upper gap as well as in the lower gap region. Then 



0=(p duj {F+{uj) + F-{uj)) 

l\u\=l 



/ du {Am5{uj -u)) + 2 / dujF{uj) + 2 duj F+{uj) (b 19) 

dio F{uj) + / dtj F{uj) 



47ri 



h+(u 



+ 2 



J^ug 



Igs 



ih+{u) 



So by multiplying ^ ^^'i^ at the last line of the above equation ( B.19 ), we can see ( |B.12| ) 



ih+{u) 
47r2A 

Finally we have shown ( |3.78| ). 

B.2.2 Proof of (|^) 

Let us prove the formula ( 3.83| ) 



Lug 



doj F{uj) + I duj F(. 



27rA 



(B.20) 









(B.21) 



Because j^j^ is holomorphic everywhere except the cuts region, from the Cauchy's 
theorem, we can see 







du} 



1 



doj 



1 



TH=00 H^) JCcuts K^) 

From T-r^ — )• — r at \lij\ — )• oo, we can also see 



^^ 1T~\ 

\u\=oo H^) 



Hence we can rewrite 



1 






1 



(B.22) 



(B.23) 



d(jj 



iccuts H^) j\u,\=i+\^\ H^) 7ic.|=i-|.r H^) 



1 



duj 



1 



+ 



|cj|=i+|e| h{uj) 7|,^|=i_|,| h{uj) 



duo 



1 



du) 



1 1 

+ 



1^1=1 V^"^(^) h (^) 



J^ug! 



K^) Jl,,, h{^) 



Iga 



(B.24) 



From first line to second line, we have used the fact that j-t-^ is holomorphic inside the 



circle IcjI = 1 — lei 



dbj 



\ui\ = l~\e\ K^) 



h(u)) 



\uj\ = l-\e\ "-l^j 



(B.25) 
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1 



In the last step from the third to the forth hne of ( B.24 ), we have used that ^4, n — — ^_, >. 
in the cuts region while 



the last hne of ( |B.24| ), 



h+{ui) 



— j^jj^ in the gaps region. Hence we can see from 
1 



h (uj) h{uj) 

1 






^ hi \ 



(B.26) 



This completes the proof of ( 3.83 ). 
B.2.3 Proof of (|^) 



Let us prove the formula (13.87] ) . 



Since j-f^ is holomorphic everywhere except the cuts, from the Cauchy's theorem, we 



h{uj) 



can see 



1 

27rl 



du 



UJ 



|a;|=oo 



h{uj) 



1 

2wi 



Because j^ 



Hence 



1 

2TTi 



7:71 



^cut 



at \uj\ — )• 00, it becomes 



duj 



UJ 



h{uj) 



dijj 



UJ 



1 

27ri 

1 



duj 



UJ 






J- r\e'K 



|l^|=oo H^) '^■^iJo W 



(B.27) 



(B.28) 



duJ -r-, — T : 



dhJ 



UJ 



Ccuts K^) 27ri 7|^|=i+|e| h{uj) 27ri J\uj\=i-\e\ H^) 

If , UJ If , UJ 

~27ri 



I'^Ni+kl 



dUJ -ri~T "I 

h[uj) 2-111 



UJ 



iri It h(uj) 



1 



^^ JUgs '^ ^(^) 



'^ hi \ 

UJ 



(B.29) 



From first line to second line we have used the fact that ^7^ is holomorphic inside the 
circle \uj\ = 1 — |e|. And from the second line to third line we have used the fact that 

[■mY = - [W)Y °^ the cuts region while (j^Y = (^)" = (^) in the gaps 
region. From the last line of ( [B.29| ), we can see 



_ 1 r UJ If UJ 



(B.30) 



This completes the proof of the formula ( |3.87D 

C. Behavior of eigenvalue distribution at large C, 

We can expect that at C — ^ co with fixed A, the eigenvalue densities would behave as 



p{a) 



1 



27rA 




\a\ < vrA) 
\a\ > vrA). 



(C.l) 



First we will verify that at C = cc with fixed finite A, the eigenvalue density in each CS 
matter theory must be this function. And after that we will elaborate the behavior of the 
eigenvalue densities at the limit more. 
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C.l Proof of the universal distribution of eigenvalue density at (^ = oo 

To verify it, it is helpful to see T(a, 6) = oo =^ a = 6, at C = oo with finite A. Here T(a, b) 
is defined in ( 0.16|) . Please note since 

^ < , \ (C.2) 



sm^ I - sin^ | cos f Jsin^ | - sin^ | 



for every |a| < a < b, it would be 

1 r 

T(a,b)<^ / da- 



sin^ I - sin2 § -^-^ cos f y'sin^ § - sin^ | 2 cos f ^sin^ | - sin^ | 

(C.3) 
From ( |C.3| ) we can see that T = oo =^ a = 6. Also in the case of a = vr, because of the 
constraint a<b<TT,a = b is automatically satisfied. So we can conclude that 

T{a,b) = oo^a = b. (C.4) 

In following subsubsections, we will verify the statement by using (|C. 



C.1.1 Regular fermion case 

First we will consider the first line of ( 3.16| ). If a positive number c is finite, y^'-'{a,b,c) 



is kept finite. Then if the positive number c is kept finite at the limit C — )• oo, A : fix, the 
RHS of ( 3.1(j ), C,y^''' {a,b,c) becomes infinite. So if positive c is kept finite at the limit, 



T(a, b) of the LHS of ( |3.16| ) also becomes infinite. Therefore if positive c is kept finite at 



the limit, it must be o = 6 at (^ = oo due to (C.4), and then the eigenvalue density becomes 
the universal distribution ( |C.1| ). 

So to complete the proof we should show that c is kept finite at the limit. Equation 
to determine the value c is 

c = fp;x{c) = — / da/)(a)log(l + 2cosae~^ + e"^^). (C.5) 

(t — Aj J_^ 

Here the function fp-x{c) is positive by the property of /O as a probability function and by 
the property that p{a') < p{a) if vr > \a'\ > \a\. At positive c we can see the following 
inequality 

2A 
0</p;A(c)<^^-^log2. (C.6) 



Solution of (C.5) is represented by the the intersecting point of the line y = x and y 



fp,\{x) in the x — y plane, which is denoted as {x,y) = {c,c). From the inequality ( C.6 ), 

we can immediately see that 

2A 
S<i^^log2^ (C.7) 



2A i„„o ;c «■„;+„ ^,,.^v,+;+-,. -f^,. „,,^v,r n ^ \ ^ 1 12 



.._,s log 2 is finite quantity for every < A < 1. So we can see that c is a finite positive 
quantity which is not infinite, then we have shown that the eigenvalue density becomes 
universal distribution (CT) at C = oo. 



^^We can not apply this analysis for the case A = 1. But in case of A = 1, the eigenvalue density is always 
p{a) = 2^, which already obeys the universal density. So we do not have to mind about the A = 1 case. 
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C.1.2 Critical boson case 



We can show in a similar way to the regular fermion case. First we will consider the first 
line of ( p.47| ). If a positive number o is finite, y^-^ia^ b, a) is kept finite. Then if a does not 
blow up to positive infinity at the limit ( — )■ oo, A : fix, the RHS of ( 3.47 ), C,y'^'''{a,b,a) 
becomes infinite. So if positive a is kept finite at the limit, T(a, 6) of the LHS of ( |3.47] ) 
also becomes infinite. Therefore if positive a is kept finite at the limit, it must he a = b 
at C = oo due to (C.4), and then the eigenvalue density becomes the universal distribution 

O- 

So to complete the proof we should show that a is kept finite at the corresponding 
limit. Equation to determine the value a is 



/7r 
da p{a) log (l — 2 cos ae~^ + e~^^) . 
-7r 



(C.8) 



Solution of (|C.j ) is represented by the the intersecting point of the line y = x and y = gp{x) 
in the x — y plane, which is denoted as (x, y) = {a, a). Please note at x > 0, 



(1-e-) 



0<5p(a;)<-21og(l-e 

Suppose that x = x > is the solution of the equation 

x = -21og(l-e-^) 

from the inequality ( |C.6| ), we can immediately see that 

. . 1 /3W5^ 
a < X = log 



(C.9) 



(CIO) 



(C.ll) 



So we can see that c is a finite positive quantity which is not infinite, then we have shown 
that the eigenvalue density becomes universal distribution ( p.l| ) at C = co- 



C.1.3 SUSY CS matter theory case 

We can show in a similar way to the previous cases. As in the previous cases, we will 
consider the first line of ( |4.16| ). If a positive number c is finite, y^^{a,b,a) is kept finite. 
Then if c does not blow up to positive infinity at the limit C — t- oo, A : fix, the RHS of 
( 4.16 ), (y'^''^{a,b,a) becomes infinite. So if positive c is kept finite at the limit, T(a, 6) of 
the LHS of (|4.16| ) also becomes infinite. Therefore if positive c is kept finite at the limit, 
it must be a = 6 at C = oo due to ( |C.4|) , and then the eigenvalue density becomes the 



universal distribution (C.l). 

So to complete the proof we should show that c is kept finite at the corresponding 
limit. Equation to determine the value c is 

1-1-9 nns n>p~'^ -A- p~'^^ 

dap{a) log 



A 



1 — 2 cos ae '^ + e 
For any c>0, 0<A<1, we can see the following inequality. 



-2c 



A 



dap{a) log 



1 + 2 cos ae '^ + e 



-2c 



1 — 2 cos ae '^ + e' 



-2c 



< 2 



log 



(C.12) 



(C.13) 
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Please note that the last term is independent of (. Our c is obtained as the solution 
X = c > of the following equation, 



X = X 



dap{a) log 



1 + 2 cos ae ^ + e 



-2x 



1 — 2 cos ae ^ + e ^^ 
Suppose that x = Cq > is the solution of the equation 

2 



log 



1 



due to the inequality ( p.l3 ), we can see 

< c< 5o = log(3 + 2\/2). 



(C.14) 



(C.15) 



(C.16) 



Hence we conclude the positive c is kept as finite quantity even if C = oo, and then we have 
shown that the eigenvalue density becomes the universal distribution ( p.l[) at C, = oo. 

C.2 Behavior of eigenvalue density at large Q, with 1 » 4 > 



We will elaborate how the eigenvalue density deviates from the universal distribution ( C.l ) 
if we gradually decrease the temperature from ^ = oo. 



C.2.1 Regular fermion theory 

To consider the behavior in the regular fermion theory, we should evaluate behavior of the 



combination of (|3.16|) , ( {3.17 ) and (|3.3|) at large Q where 6 ~ a + e with e <C 1. (As in 
section C.2 in [y.) At this limit, these three equations are evaluated asymptotically as 



1 



1 



27rA V sin(a) 



log 



8sin(a) 



e 
V 



cos(a) 



sin^fa' 



+ e 



1 + log 
ysva.a 



8sin(a) 



+ 0\e' 



cosh y + cos a 2(coshy + cosa) 



(C.17) 



and 



1 / / \ 1 / »sm(a) 

2 a + 2 cot (a) log 



27rA 



l + ^Tdy r'"" . 

2tt J^ cosh y + cos a 4tt 



sin (a) 



-1 + log 



8sin(a) 



+ 0[e'] 



C f°^ 1 y cosh y sin a 
e^- / dy 



(CIS) 



(cosh 2/ + cos a) 2 ' 



and 



1 



2^(1 - A) J_^^ 
respectively. From these we can obtain e and a, b as 



dalog(l + 2cosae " + 6 ) 



a • / \^ / sin(7rA) 
e = 8sm(7rA)exp ( AC 



dy 



y 



cosh y + cos vrA 



+ 



a = ttA- -e. 



vrA + i. 



(C.19) 

(C.20) 
(C.21) 
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Since 

Jc cosh y + cos vr A 

C is estimated as C ~ — log e from (|C.20D . 

Let us consider the asymptotic behavior of the eigenvalue density ( 3.18| ) in e — )• limit. 
Here p2,tg of ( p. 18 ) is the same functional form as the one in the GWW type model (7.6) of 
||l|, while /9^{ of ( 3.1^ ) is regarded as an additional term appearing in the regular fermion 



theory. We can see that p^'{ falls off as p^L ~ — eloge — )■ 0. On the other hand p2,tg is 
kept finite. So at the limit, p becomes 

11 1 

p(a) = P2.tg{>^, ttA - -e, vrA + -e; a) = -^— cos"^ ^/oT (C.23) 

z z vr^A 

where a = a + aie where < ai < 1. This is very similar behavior to the one in the GWW 
type model, (7.11) of Q. 

Although the functional form of p at large ( is almost same as the one in the GWW 
type model, the range of e in the regular fermion theory would be different from e in the 
GWW type model. To evaluate the e, we should calculate the factor Yrf in ( p. 20 ), and 



compare with e = 8 sin(7rA) exp ( —^—j-^X(^ 1 in the GWW case. By numerical calculations, 
we can confirm that 

Yrf > 1 (C.24) 

for every value of A. (You can see the plots of this quantity Yrf in Fig. 11(a)). Hence we 



can see that e at (A, () of the regular fermion theory (C.20) is smaller than e at the same 



(A, () of GWW model. This means that the shape of the eigenvalue density function of the 
regular fermion theory has a sharper cliff than the one in GWW model. (Schematic graph 
of the eigenvalue density is depicted at Fig. 10(a)| .) 



C.2.2 Critical boson theory 

Also in the critical boson theory, we consider how the eigenvalue density deviates from 
universal configuration ( p.l| ) as we decrease the temperature. 

To consider the behavior, we should investigate the combination of the conditions 
( 3.47 ), (|3.48| ), and ( ^3.36 ) in the large C where b = a + e with e ^ 1. At this hmit, the 



asymptotic behavior of these three equations are 

1 f^ 1 ^^j8M^\_^co^f^^^^j8M^^^^^^., 



27rA \ sin(a) \ ^ J sin (a) 
C /"°° , f y ysina 



dy 



2-71 J ^ V *^osh y — cos a 2 (cosh y — cos a) 

and 



(C.25) 



1 / , ,, /8sin(a)\ 1 / , /8sin(a) , ■ ■ ^. ., 
' 2 a + 2 cot(a) log ( -^^ ) - e—^-j { -1 + log ( j^^ ) ) ) + 0[e' 



27rA \ V f / sin^(a) 

C /"°° r y cos a (" f°° y cosh y sin a 



27r Jo- cosh y — cos a 4tt J^ (cosh y — cos a)^ ' 



(C.26) 
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(a) 

Figure 10: Plots of cos^^(y/ar) curves, which are the eigenvalue density functions on the cut 
domain ttA— |<Q;<7rA+|atC— > cx). (We did not plot precise e here.) By using these schematic 
graphs, we are trying to compare the behaviors of the eigenvalue densities; the one of the GWW, 
the one of the regular fermion (the critical boson) and the one in the SUSY CS matter theory. 
Dotted blue line is the one in the GWW model and the bold line is the regular fermion's one, and 
the dashed red line is the one in the SUSY CS matter theory. The regular fermion theory has a 
sharper slope than the one in the GWW model, moreover the SUSY theory has a sharper slope 
than the one in the regular fermion theory. At fixed A, the center of the cut domain is always fixed 
at A. 



and 



1 



/ttA 
da log (1-2 cos ae~" + e''^") , 
-ttA 



2vrA j_^;^ 
respectively. From these, e, a and b are evaluated as 

/I f° 

e =8sin(7rA) exp I — -sin(7rA)CA / 



dy y 



cosh y — cos vrA 



a =7rA 



Since 



K 



cb 



6 = 7rA + 



dy 



y 



cosh y — cos TT A 



0{e% 



(C.27) 

(C.28) 
(C.29) 

(C.30) 



Q is evaluated as C ~ — loge from (C. 



Let us consider the asymptotic behavior of the eigenvalue density ( |3.49| ) in e — )• limit. 
Here p2,tg of ( p.49| ) is the same functional form as the one in the GWW type model (7.6) 
of PI, while p'^\ of ( |3.49 ) is regarded as an additional term appearing in the critical boson 
theory. We can see that p^'j falls off as p^^^ ~ — eloge — )■ 0. On the other hand p2,tg is 
kept finite. So at the limit, p becomes 

111 

p{a) = P2.tg{\ ttA - -e, vrA + -e; a) = ^^ cos"^ y/a{. (C.31) 



7r2A 



This is also very similar behavior to the one in the regular fermion case as well as the one 
in the GWW type (7.11) of §. 

The range of e in the critical boson theory would be also different from e in the GWW 
type model. To evaluate the e, we should calculate the factor Yd, in ( C.28| ), and compare 
with e in the GWW case. By numerical calculations, we can confirm that 



>;fe>i 



(C.32) 
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for every value of A. (See the plots of Y^b in Fig. 11(a) ). This means that the shape of the 
eigenvalue density function of the critical boson theory has a sharper cliff than the one in 
GWW model. (Schematic graph of the eigenvalue density is depicted at Fig. 10(a)| .) We 
can see that e at ( p.28| ) as well as Y^b are dual to e at ( C.20D and Y^j under the level-rank 
duality. We can see the duality between Y^b and Y^f also from the graph Fig. 11(a) . 



C.2.3 Supersymmetric CS matter theory 

Similar to the other theories cases, we should consider the behavior of the combination of 



( 4.16 ), ( 4.17 ) and (4^) at large ("• At the limit these behave as 



1 



27rA 



1 



smla 



■log 



8sin(a) 






2y cosh y 



and 



27rA 



cos(a) / ^ 
■e-^ -1 + log 
sm^(a) V 

2y sin a cos a cosh y 

\ cosh^ y — cos^ a (cosh^ y — cos^ a)^ 

8sin(a) 



ism a 



+ 0[e2 



(C.33) 



2 a + 2 cot (a) log 



1 



sin (a) 



1 + log 



8sin(a) 



+ 0[e' 



C f°° , 2w cos a cosh y 

=1 + 7^ / rfy 



2lT 



cosh y — cos^ a 47r 



e( f°° 2y coshy sina(cosh^ y + cos^ a) 
IT dy — 



(cosh y — cos^ a) 



and 



1 

2^ 



ttA 



da log 



1 + 2 cos ae ^ + e 



-2c 



-ttA 



1 — 2 cos ae '^ + e 



-2c 



From these, e, a and h are evaluated as 



/ , X / sin(7rA) ^, 
8sm(7rA)exp ( tt^CA 



dy 2ycosh.y 
cosh y — cos^ ttA 



+ 



vrA-|, 



.A + -. 



Since 



Y 



dy 2y cosh y 
cosh y — cos^ ttA 



0{e'), 



(C.34) 
(C.35) 

(C.36) 
(C.37) 

(C.38) 



we can see (^ ~ — log e. 

Let us consider the asymptotic behavior of the eigenvalue density ( 4.15| ) in the limit. 
Here p2,tg of ( [4.15| ) is the same functional form as the one in the GWW type model (7.6) 
of Q, while pI^^J' of ( 4.15 ) is regarded as an additional term appearing in the SUSY CS 
matter theory. We can see that Pi"!^ falls off as Pi^^^ ~ — eloge — ?• 0. On the other hand 
P2,tg is kept finite. So at the limit, p becomes 



p{a) = /02.tg(A,7rA 



1 , 1 ^ 1 

-e,7rA + -e;a) = ^-r cos 

2 2 vr^A 



ai. 



(C.39) 



This is almost same behavior as the other CS matter theories. 
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Figure 11: The points with black filled circle are plots oiY = Ygu in the SUSY CS matter theory. 



The circle and square plots without being filled are Y 
fermion theory and Y = Yd, — J dy 



K 



rf 



J c ^ c 



. -^ r in the regular 

cosny+cosTrA ^ 

at the critical boson theory respectively. Dotted 



cosh y— cos ttA 

horizontal line indicates Y = 1. We can see that for any A, Ygu > Y^f > 1 and K;„ > Ycb > 1 



The range of e in the SUSY CS matter theory would be the smallest among the 
CS matter theories which we have already studied in this paper. We can estimate e by 
calculating the factor Ysu in ( p.36| ). (The plots of Ysu is in Fig. 11(a) ). From this, we 
can see Ysu > Ycb-, (Xrf) > 1- This means that e in the SUSY CS matter theory is the 
smallest and shape of the eigenvalue density function of the SUSY CS matter theory has 
the sharpest cliff among the CS matter theories. (Schematic graph of the eigenvalue density 
is depicted at Fig. 10(a) .) We can see that e at ( p.36[ ) as well as Ygu are self-dual under 



the level-rank duality. We can see the self-duality of Ygu also from the graph Fig. 11(a) 



References 

[1] S. Jain, S. Minwalla, T. Sharma, T. Takimi, S. R. Wadia, et. 



Phases of large N vector 



Chern-Simons theories on S'^ x S^ , 



arXiv: 1301. 6169 



[2] I. Klebanov and A. Polyakov, AdS dual of the critical 0(N) vector model, Phys.Lett. B550 



(2002) 213-219, |hep-th/0210114 



[3] M. A. Vasiliev, Consistent equation for interacting gauge fields of all spins in 
(3+l)-dimensions, Phys.Lett. B243 (1990) 378-382. 

[4] M. Vasiliev, Nonlinear equations for symmetric massless higher spin fields in (A)dS(d), 
Phys.Lett. B567 (2003) 139-151, |iep-th/0304049 |. 



[5] S. Giombi and X. Yin, The Higher Spin/Vector Model Duality, arXiv: 1208.4036 



[6] S. Giombi, S. Minwalla, S. Prakash, S. P. Trivedi, S. R. Wadia, et. ai, Chern-Simons Theory 



with Vector Fermion Matter, Eur.Phys.J. C72 (2012) 2112, [irXiv: 1110 .4386|. 



[7] C.-M. Chang, S. Minwalla, T. Sharma, and X. Yin, ABJ Triality: from Higher Spin Fields to 
Strings, [arXiv: 1207 .4485 . 



59 



E. Witten, Anti-de Sitter space, thermal phase transition, and confinement in gauge theories, 



[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 



Adv.Theor.Math.Phys. 2 (1998) 505-532, |lhep-th/9803131 



O. Aharony, S. Giombi, G. Gur-Ari, J. Maldacena, and R. Yacoby, The Thermal Free Energy 
in Large N Chern-Simons-Matter Theories, 



arXiv: 1211.4843 



M. R. Douglas and V. A. Kazakov, Large N phase transition in continuum QCD in 



two-dimensions, Phys.Lett. B319 (1993) 219-230, [liep-th/9305047 



X. Arsiwalla, R. Boels, M. Marino, and A. Sinkovics, Phase transitions in q-deformed 2-D 



Yang-Mills theory and topological strings, Phys.Rev. D73 (2006) 026005, |hep-th/0509002 |. 



N. Caporaso, M. Cirafici, L. Griguolo, S. Pasquctti, D. Seminara, et. al.. Topological strings 
and large N phase transitions. L Nonchiral expansion of q-deformed Yang-Mills theory, JHEP 



0601 (2006) 035, | |hep-th/0509041 |. 

D. JafFeris and J. Marsano, A DK phase transition in q-deformed Yang-Mills on S**2 and 



topological strings, hep-th/0509004. 

A. Kapustin, B. Willett, and I. Yaakov, Nonperturhative Tests of Three- Dimensional 



Dualities, JHEP 1010 (2010) 013, arXiv: 1003.5694 



A. Kapustin, B. Willett, and I. Yaakov, Tests of Seiberg-like Duality in Three Dimensions, 



arXiv: 1012.4021 



O. Aharony, G. Gur-Ari, and R. Yacoby, d=3 Bosonic Vector Models Coupled to 
Chern-Simons Gauge Theories, JHEP 1203 (2012) 037, 



arXiv: 1110.4382 



O. Aharony, G. Gur-Ari, and R. Yacoby, Correlation Functions of Large N 



Chern-Simons-Matter Theories and Bosonization in Three Dimensions, arXiv: 1207.4593. 



E. Ivanov, Chern-Simons matter systems with manifest N=2 super symmetry, Phys.Lett. 
B268 (1991) 203-208. 

D. Gaiotto and X. Yin, Notes on superconformal Chern-Simons-Matter theories, JHEP 0708 



(2007) 056, |arXiv:0704.3740 



R. D. Pisarski and S. Rao, Topologically Massive Chromodynamics in the Perturbative 
Regime, Phys.Rev. D32 (1985) 2081. 

W. Chen, G. W. SemenofF, and Y.-S. Wu, Two loop analysis of nonAbelian Chern-Simons 



theory, Phys.Rev. D46 (1992) 5521-5539, |hep-th/9209005 |. 

J. Maldacena and A. Zhiboedov, Constraining Conformal Field Theories with A Higher Spin 



Symmetry, arXiv: 1112 . 1016. 

J. Maldacena and A. Zhiboedov, Constraining conformal field theories with a slightly broken 
higher spin symmetry. 



arXiv: 1204.3882 



S. Banerjee, S. Hellerman, J. Maltz, and S. H. Shenker, Light States in Chern-Simons Theory 
Coupled to Fundamental Matter, 



arXiv: 1207.4195 



S. Jain, S. P. Trivedi, S. R. Wadia, and S. Yokoyama, Supersymmetric Chern-Simons 
Theories with Vector Matter, JHEP 1210 (2012) 194, farXiv: 1207.475C| ]. 



S. Yokoyama, Chern- Simons- Fermion Vector Model with Chemical Potential, 
arXiv:1210.410g| . 



60 



[27; 

[28 
[29 

[3o; 

[31 

[32 
[33 

[34; 

[35 

[36; 

[37; 

[38 
[39; 

[4o; 

[41 

[42 
[43 



S. Banerjee, A. Castro, S. Hellerman, E. Hijano, A. Lepage- Jutier, et. al., Smoothed 
Transitions in Higher Spin AdS Gravity^ arXiv : 1209 . 5396| . 



G. Gur-Ari and R. Yacoby, Correlators of Large N Fermionic Chern-Simons Vector Models^ 
arXiv: 1211 .18661 . 



S. H. Shenker and X. Yin, Vector Models in the Singlet Sector at Finite Temperature, 
arXiv:1109.35ig| . 



M. Blau and G. Thompson, Derivation of the Verlinde formula from Chern-Simons theory 



and the G/G model, Nucl.Phys. B408 (1993) 345-390, [ tiep-th/9305010 1 



D. Gross and E. Witten, Possible Third Order Phase Transition in the Large N Lattice Gauge 
Theory, Phys.Rev. D21 (1980) 446-453. 



S. R. Wadia, A Study of U(N) Lattice Gauge Theory in 2-dimensions, arXiv: 1212.2906. 



S. R. Wadia, N — infinity phase transition in a class of exactly soluble model lattice gauge 
theories, Phys.Lett. B93 (1980) 403. 

B. Sundborg, The Hagedorn transition, deconfinement and N—4 SYM theory, Nucl.Phys. 



B573 (2000) 349-363, [iep-th/9908001 



O. Aharony, J. Marsano, S. Minwalla, K. Papadodimas, and M. Van Raamsdonk, The 
Hagedorn - deconfinement phase transition in weakly coupled large N gauge theories. 



Adv.Theor.Math.Phys. 8 (2004) 603-696, |lhep-th/0310285 | 



O. Aharony, J. Marsano, S. Minwaha, K. Papadodimas, and M. Van Raamsdonk, A First 
order deconfinement transition in large N Yang-Mills theory on a small S**3, Phys.Rev. D71 



(2005) 125018, [ lhep-th/0502149 |. 

O. Aharony, J. Marsano, and M. Van Raamsdonk, Two loop partition function for large N 



pure Yang-Mills theory on a small S**3, Phys.Rev. T>74 (2006) 105012, hep-th/0608156 



K. Papadodimas, H.-H. Shieh, and M. Van Raamsdonk, A Second order deconfinement 
transition for large N 2+1 dimensional Yang-Mills theory on a small two-sphere, JHEP 0704 



(2007) 069, [ |hep-th/0612066 |. 

M. Mussel and R. Yacoby, The 2-loop partition function of large N gauge theories with 



adjoint matter on S**3, JHEP 0912 (2009) 005, [jarXiv : 0909 . 0407[| . 

L. Alvarez-Gaume, C. Gomez, H. Liu, and S. Wadia, Finite temperature effective action, 



AdS(5) black holes, and 1/N expansion, Phys.Rev. D71 (2005) 124023, [lhep-th/0502227 |. 



L. Alvarez-Gaume, P. Basu, M. Marino, and S. R. Wadia, Blackhole/ String Transition for the 
Small Schwarzschild Blackhole of AdS(5)x S**5 and Critical Unitary Matrix Models, 



Eur.Phys.J. C48 (2006) 647-665, [ tiep-th/0605041 1 



M. R. Douglas, Chern-Simons- Witten theory as a topological Fermi liquid, liep-th/9403119 



A. Giveon and D. Kutasov, Seiberg Duality in Chern-Simons Theory, Nucl.Phys. B812 
(2009) 1-11, |arXiv:0808.0360|| . 



[44] F. Benini, C. Closset, and S. Cremonesi, Comments on 3d Seiberg-like dualities, JHEP 1110 



(2011) 075, [|arXiv: 1108. 5373(1 



61 



